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ABSTRACT 

Intergalactic medium (IGM) metal absorption lines observed in z > 6 spectra offer 
the opportunity to probe early feedback processes, the nature of enriching sources, and 
the topology of reionization. We run high-resolution cosmological simulations including 
galactic outflows to study the observability and physical properties of 5 ions (C II, C iv, 
Oi, Sin, Siiv) in absorption between z = 8 — » 5. We apply three cases for ionization 
conditions: Fully neutral, fully rcionizcd, and a patchy model based on the flux from 
the nearest galaxy. We find that our simulations can broadly fit available z ~ 5 — 6 
IGM metal-line data, although all observations cannot be accommodated with a single 
ionization condition. Variations in O I absorbers among sight lines seen by Becker et al. 
(2006) suggest significant neutral IGM patches down to z ~ 6. Strong C iv absorbers 
at z ~ 6 may be the result of ionization by the galaxy responsible for that enrichment, 
although the identification of the neighboring galaxy will have to wait to confirm 
this. Our outflows have typical speeds of ~ 200 km/s and mass loading factors of 
~ 6. Such high mass loading is critical for enriching the IGM to the observed levels 
while sufficiently curtailing early star formation to match the observed rest-frame UV 
luminosity function. The volume filling factor of metals increases during this epoch, 
but only reaches ~ 1% for Z > 10~ 3 Z Q by z — 5. Detectable absorbers generally trace 
inhomogeneously-distributed metals residing outside of galactic halos. C iv is an ideal 
tracer of IGM metals at z ~ 5 — 6, with dropping global ionization fractions to either 
higher or lower redshifts. This results in a strongly increasing global C iv mass density 
(fi(Crv)) from z = 8 — * 5, in contrast to its relative constancy from z = 5 — > 2. Our 
simulations do not support widespread early IGM enrichment from e.g. Population III 
stars, as this would overpredict the numbers of weak C iv absorbers in the latest data. 
High-z absorbers arise from metals mostly on their first outward journey, at distances 
5-50 physical kpc, and often exhibit broad profiles (Sv > 200 kms -1 ) as a result of 
outflowing peculiar velocities in the strongest systems. The galaxies responsible for 
early IGM enrichment have typical stellar masses of if) 70-8 - 5 M , and star formation 
rates < 1 M /yr. Future facilities will be able to study the high-z galaxy-absorber 
connection in detail, revealing a wealth of information about feedback processes in the 
reionization epoch. 

Key words: intergalactic medium, galaxies: formation, galaxies: high-redshift, early 
Universe, cosmology: theory, methods: numerical 



1 INTRODUCTION 

The Universe undergoes its last major transition at z ~ 
6 as the first stars and early galaxies finish the process 
of reionizing the previously neutral intergalactic medium 
(IGM). Observations have only recently cracked the < 1 
Gyr Universe w ith the discovery o f galaxies tracing early 



Uyr Universe w ith the discovery o l galaxies tracing early 
star formation dBunker et al.l 120041: [D ickinson et al.l 2004; 
lYan fe Windhorsdl2004l : iBouwens etali 2006), quasars pro- 



viding ionizing photons (|Fan et al.ll200ll , l2003l ). and the com- 
plete absorption b y the I GM of all Lyman-a photons, i.e. the 
iGunn fc Peterson! (|l965h trough, in z > 6 quasar spectra 
(iFan et alj|2002f). 



A new type of high-z observation is the recent detection 
of metal-line absorbers in the most distant quasar spectra, 
probi ng the nucleosynth etic products of early star forma- 
tion . [Becker^etjiL (2006, hereafter BSRS) surveyed 9 z > 5.5 
quasars and found 4 O I systems in the most distant ob- 
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ject in their sa mp le, J1148+5251. iRvan- Weber et all |2006), 
ISimcoel i|2006l ). fc lRvan- Weber et al l ([20091 ) discovered mul- 
tiple Civ absorbers at similar redshifts (2 ~ 6). Despite the 
metals presumably being formed within galaxies, these ab- 
sorbers are most readily explained as arising in the diffuse 
IGM. But such absorbers are rare and perhap s highly inho- 
mogeneous: For instance, iBecker et al.l (|2009l ) finds no C IV 
along four lines of sight (LOSs) observed at higher signal- 
to-noise {S/N) between z — 5.3 — 6.0. 

The simple fact that metals exist in the early IGM begs 
the question, how did they get there? Population III (i.e. 
metal-free) and early stars have been proposed to enrich a 
significant volume of the early Universe, as early haloes har- 
boring primordial galaxies have small physica l scales and 
shallow potential wells (e .g. iMadau et al.l [200ll : iGreif et alj 
120071 ; IWise fc Abell 12008). The relat ively invariant amount 
of C iv absorptio n between z ~ 2 — 5 |Songailall200ll) an d ex- 
tended to z ~ 6 l|Rvan- Weber et al.ll2006l ; lsimcoell20of3 ) has 
been interpreted to mean that t he majority of IGM met - 
als were injected at z > 6 (e.g. IScannapieco et aL 2002). 
But the most recen t data set from IBecker et al.l ( 20091 ) and 
I Ryan- Weber et all (|2009l ) shows a decline of a factor of 
~ 4x at z > 5.3 compared to z < 4.7, suggesting that 
IGM metallicity is increasing at early epochs, assuming 
typical ionization conditions do not evolve much. Further- 
more, z ~ 3 Civ absorbers are spatially correlated with 
Lyma n Break Galaxies (LBG s ), sug gesting ongoing enrich- 
ment l|Adelberger et alj l2003t 120051 ). although these mea- 
surements are also consistent with metals being injected by 
dwarf galaxie s at z = 6— 12 coupled wi th subsequent growth 
of clustering l|Porciani fc Madaull2005l ). 

lOppenheimer fc Pavel l|2006l . hereafter OD06) explored 
C IV absorbers in cosmological hydrodynamic simulations, 
employing various heuristic prescriptions for galactic out- 
flows to enrich the z ~ 2 — 6 IGM. They determined that pre- 
scriptions based on momentum-driven winds were best able 
to reproduce Civ observations. OD06 also showed that an 
increasing ionization correction for C IV, from both cosmic 
evolution and energy injection from outflows, could mask a 
comparable increase of ~ x 10 in the true IGM metallicity 
from z ~ 6 — > 2. Hence widespread early enrichment from 
some exotic stellar population is not required to explain the 
relatively invariant C IV absorption over this redshift range. 

Another important question is, how do early galax- 
ies relate to the IGM metal-l i ne abs orbers at z > 6? 
iDave. Finlator. fc Qppenheimerl l|2006l . hereafter DFO06) 
explored reionization-epoch galaxies in the simulations of 
OD06, finding that strong outflows are required to match 
the obser ved luminosity funct io n (LF) of z ~ 6 i- band 
dropouts (|Bouwens et al.ll2006l ). iFinlator et~aH l|2007h fur- 
ther found that the star formation histories from outflow 
simulations can reproduce broad-band spectral energy dis- 
tributions (SEDs) of gala xies such as the le nsed z ~ 6.7 
galaxy Abell 2218 KESR l|Kneib et al.ll2004l ). The success 
of these models lies in the high mass loading of early 
galactic superwinds, which pushes the majority of high- 
z metals nucleosynthesized i n massive stars into the IGM 
|Dave fc Qppenheimerl 120071 . Oppenheimer & Dave 2008, 
hereafter OD08). Typically, many times more mass is driven 
out of galaxies than forms into stars. Hence if our models 
are correct, observations of IGM metal-line absorbers may 



be the best way to account for the majority of metals in the 
z > 6 Universe. 

Finally, what can IGM metal-line absorbers reveal 
about the topology of reionization? Fifth year data 
from the Wilkinson Microwave Anisotropy Probe (WMAP, 
iHinshaw et al.ll2008h indicate the Universe was reionized at 
z ~ 11 ± 1.4 if it proceeded instantaneously. Numerical sim- 
ulations suggest that reionization is a protracted process 
whereby individual sources create Hn bubbles at z > 10 
that propagate outward until overlap by z ^ 6 (e.g. | Gnedinl 
200d; iFurlanetto fc Ohl 120051 ; lllievl I2OO6I ; iLidz et al l 12007 



Finlator et aLl F2009 ) . Observations of H I absorption are di - 
vided over when reionization ends, with iFan et al l (|2006ri 
showing a discontinuity in the IGM ionization rate at z ~ 6 
suggesting individ ual ionization bubble s have finally finished 
overlapping, while IBecker et al.l (|2007l ) finds that the Lya 
optical depth distribution shows no evidence of a sudden 
end to reionization at this epoch. Metal lines provide an ad- 
ditional handle on how reionization progresses, allowing the 
bulk of cosmic volume to be probed in absorption despite a 
fully saturated Gunn-Peterson trough. 

In this paper we consider the above questions by 
studying z = 5 — 8 metal-line absorbers in cosmologi- 
cal hydrodynamic simulations that incorporate our favored 
momentum-driven wind model. In this wind model, the 
outflow speed and mass loss rate are tied to galaxy mass 
using a heurist ic prescription as expected for momentum- 
driven winds l|Murrav. Quatert. fc Thompson! 12005). This 
wind model has enjoyed a number of successes when com- 
pared to a wide range of observations at z < 6; most 
relevant is that it fairly uniquely matches observations of 
z ~ 2 — 5 Civ absorbers (OD06). The scaling of wind 
speed with galaxy velocity in this model fo llows the rela- 
tion observed in lo cal starburst outflows (e.g. Martin 2005a; 
Rupke et al]|2005l ). We note that such scalings may arise 
from other physical scenari os besides momentum (or rad i- 
ation) driven outflows (e.g. iDalla Vecchia fc Schavel 120081 ) : 
here we merely employ the scalings, without reference to 
the underlying physical driving mechanism. 

We employ state-of-the-art Gadget-2 cosmological hy- 
drodynamic simulations with a quarter billion particles run 
to z = 5; we describe details in §2. The salient predicted 
galaxy and IGM properties are examined in §3, and com- 
pared to our findings in DFO06. In §4 we consider vari- 
ous cases for the ionization background, including the ex- 
treme cases of fully reionized and fully neutral, as well as a 
background with a variable intensity intended to represent 
ionized bubbles around individual star forming sources. §5 
presents a plethora of simulated observations which are com- 
pared to available z > 5 metal-line observations, as well as 
predictions for comparison with future data. We then study 
the physical and environmental parameters of simulated ab- 
sorbers in §6, with a focus on how absorbers relate to galax- 
ies and when they were injecte d into the IGM . §7 p rovides 
a summary. Throughout we use lAsplund et aL I (|2005l ) abun- 
dances when calibrating metallicities relative to solar. 



2 SIMULATIONS 

We employ our modified version of the N- body + Smoot hed 
Particle Hydrodynamics code Gadget- 2 (|Springelll2005h to 
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run three high-resolution 2 x 512 3 particle simulation to 
explore the high redshift Universe. Our simulations self- 
consistently enrich the IGM via galactic outflows. A com- 
plete description of our code, including our wind model, can 
be found in §2 of OD06 with further modifications described 
in §2 of OD08; however we also provide some of the fun- 
damental details involving star formation, enrichment, and 
feedback here. 



2.1 Star Formation & Chemical Production 

Star formation (SF) is modeled using a subgrid recipe where 
each gas particle above a critical density is treated as a set 
of cold clouds embedded in a warm ionized medium, sim- 
ilar to the interstella r mediu m (ISM) of our own Galaxy 
ISpringel fc HernauistJ |2003a). Gas particles eligible for SF 
undergo instantaneous self-enrichment from Type II super- 
novae (SNe). The instantaneous recycling approximation of 
Type II SNe energy to the warm ISM phase self-regulates 
SF resulting in convergence in star formation rates (SFRs) 
when looking at higher resolutions. The SFR is scaled to fit 
th e disk-surfa c e den sity- SFR observed in the local Universe 
bv lKennicuttl (|l998l fl 

Star formation below 10 Mq, assumed not to result in 
Type II SNe, is decoupled from their high mass counterparts 
using a Monte Carlo algorithm that spawns star particles. 
The metallicity of a star particle remains fixed once formed; 
however, since Type II SNe enrichment is continuous while 
stars are formed stochastically, every star particle invariably 
has a non-zero metallicity. 



2.2 Chemical Yields 

The chemical yields at high redshift are dominated by Type 
II S Ne, for which we use the metallicity-dependent yields 
from lChieffi fc Limongil 112004 ) as described in OD08. These 
yields vary less than 15% for the 4 species we track (carbon, 
oxygen, silicon, & iron) as long as Z > 10 -6 , which applies 
to effectively all SF in our cosmological simulation. We do 
not consider exotic enrichment yields fro m zero metallicity, 
Pop III, or very massive s tars (VMSs, e.g. iHeger fc Wooslevl 
2002, 2008), and assume a lChabrierl (120031 ) initial mass func- 
tion (IMF) throughout. Although it may be that such ex- 
otic stars have an effe ct on the IGM particularly in terms of 
metal-line ratios (e.g. lAguirre et al. I l2004l2008| ). the yields 
used here provide a standard baseline for comparison. The 
baseline yields are [O/C1=0 . 27 an d [Si/O]=-0.16 if we con- 
sid er thelChieffi fc Limongil l|2004 ) Z = 10~ 3 yields assum- 
ing lAselun^etjiL I |2005l ) abundances; the carbon mass yield 
ejected by SNe is 3.0 x 10 -3 the total stellar mass in a 
Chabrier IMF. The integrated gas-phase abundance ratios 
of the z = 6 dl6n256vzw si mulation box are [O / C1=0 .23 and 
[Si/O]=-0.13 reflecting the IChieffi fc Limongil i|2004h yields 
inputted in our simulations. 

To allow the reader to scale our results to a differ- 
ent set of theoretical Type II SNe yields we provide Table 
[1] listing the integrated yields assuming a IChabrierl (|2003l ) 
IMF slope with f3 = 0.18, where (5 is the stellar IMF mass 



fraction in the range indicated in the second column^. We 
have interpolated the published yields to Z = 10 -3 (i.e. 
0.08 Zq) using their various metallicity models, because this 
is near the median stellar metallicity at z = 6 (Z s t e u a r{z = 
6) = 1.48 x 10~ 3 ). There exists little spread for the oxy - 
gen and silicon yie l ds be tween IChieffi fc Limongil (l2004h , 
IWooslev fc Weaver! (|l995l ), and IPortinari et all (|l998f ) de- 
spite different SNe m ass progenitor range s . 38% less carbon 
is produced from the lWooslev fc Weaver! l| 19951 ) yield mod- 
els, which is significant when we consider a measurement 
such as fi(Clv) ( §E2M - 

We note that there may be reasons to favor a more 
top-h eavy or bottom-light IMF at high redshifts (e.g. iDavel 
2008, i.e. /3 increases). However, since our simulations are 
well-matched to observed z ~ 6 rest-UV LFs (as we will 
show in §3), and since metals mostly come from the same 
massive stars responsible for UV emission, our total metal 
production should be reasonably well-constrained. Along 
these lines we show t he yields using the 30 — 120 Mq from 
IPortinari et al.l (| 19981 ) in Table [T] silicon appears to be sig- 
nificantly under-produ ced in SNe ejecta . For comparison 
we include yields from IHeger fc Wooslevl l|2008l ) from zero 
metallicity stars. These y ields show less diffe rence relative to 
IChieffi fc Limongil l|2004h than the top-heavv lPortinari et all 



(1998) (except for oxygen), perhaps indicating just how hard 
it is to clearly identify the nucleosynthetic signatures of Pop 
III stars. 



2.3 Feedback Model 



The im plementation of winds follows ISpringel fc Hernquistl 
j|2003al ), using a Monte Carlo ejection mechanism gov- 
erned by two parameters, the outflow velocity w w i n d and 
the mass loading factor r), where the mass loading factor 
is defined as the outflow mass rate (M w i n d) in units of 
the galaxy's SFR, i.e. M w i„d — r)X SFR. The wind model 
we employ here follows scalings of « w ind and rj expected 
for momentum-driven winds; this fairly uniquely matches a 
broad range of observations including IGM enrichment be- 
tween z = 6 — > 1.5 traced by Civ (OD06) a nd at z < 0.5 
traced by Ovi (lOppenheimer fc Dave! 120091 ). the LFs of 
2 = 6 galaxies (DFO06), the galaxy mass-metallicity re- 
lations l|Finlator fc Davell2008l ) . and the enrichm ent and en- 
tropy levels in intragroup gas l|Dave et al.l [20 081 . The scal- 
ings depend on galaxy velocity dispersion (a): 



«wind = 3(7 — 1, 

3 

a 



'I 



(1) 
(2) 



where ft is the luminosity factor, which is the luminosity 
of the galaxy in units of the critical (or sometimes called 
Eddington) luminosity of the galaxy, and co provides a nor- 
malization for the mass loading factor. Here we randomly 
select /l = [1.05, 2] for each SPH particle we kick, foll owing 
observations of local starbursts bv lRupke et~all l|2005l ). and 
we take <to = 150 kms" 1 . a is estimated from the galaxy's 
mass following IMo et all (Il998l ). where the mass is obtained 
from an on-the-fly friends-of- friends galaxy finder. 



1 Even if this scaling is not applicable at z > 5, the accretion 
rate determines the SF as we argue in £|3.2I 



2 p = 0.18 for a Chabrier IMF for stars 10 - 100 Mq, which wc 
assume are the stars that undergo Type II SNe. 
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Table 1. Type II Supernovae Yields" 



Reference 


Mass Range 


Z 




[C/H] 


[O/H] 


[Si/H] 


[O/C] 


[Si/O] 


Chieffl & Limongi (2004) ft 


13-35 Mq 


io- 


3 


0.13 


0.40 


0.24 


0.27 


-0.16 


Woosley & Weaver (1995) 


12-40 Mq 


io- 


3 


-0.08 


0.39 


0.25 


0.47 


-0.14 


Portinari et al. (1998) 


10-120 Mq 


10" 


3 


-0.05 


0.43 


0.17 


0.48 


-0.26 


Portinari et al. (1998) 


30-120 Mq 


10" 


3 


-0.21 


0.35 


-1.54 


0.56 


-1.89 


Heger & Woosley (2008) 


10-100 Mq 







0.02 


0.20 


-0.29 


0.18 


-0.49 


dl6n256vzw 2 = 6 SPH 














0.23 


-0.13 



a Assuming a Chabrier IMF slope and /3 = 0.18 (i.e. 18% of stellar IMF is covered by the mass range in column 2 and only this mass 

range contributes to Type II SNe yields.). 
b We apply this set of metallicity-dependent yields in our simulations. 



2.4 Runs 

The simulations adopt the c osmological paramete rs based 
on the 5-year WMAP results (|Hinshaw et al.ll2008h . namely 
Q = 0.25, n A = 0.75, n b = 0.044, H = 70 km s _1 Mpc^ 1 , 
erg = 0.83, and n = 0.95. The value of as slightly exceeds 
the favored WMAP data alone, but agrees better with the 
combined results including Type la SNe and baryonic acous- 
tic oscillation data. A greater as results in appreciably more 
SF at high- z, which is important since SF drives the enrich- 
ing outflows. Our runs with this cosmology are named the 
d-series; the primary run is named dl6n512vzw using our 
standard naming convention, indicating the box spans 16 
comoving /i -1 Mpcs with 512 3 each of gas and dark matter 
particles. The vzw suffix indicates our implementation of the 
momentum-driven wind model as described in OD08 (with 
one inconsequential exception: we no longer impose an upper 
limit due to SN energy limitations on « w i n d). The gas par- 
ticle mass is 5.4 x 10 5 Mq, which translates to a minimum 
galaxy baryonic mass of 1.7 x 10 7 Mq assuming a galaxy is 
resolved with 32 particles (which we will demonstrate later). 
The softening length is set to 0.6 comoving ft _1 kpc. 

Two more simulations, d8n512vzw and a d32n512vzw, 
are also run to explore numerical resolution convergence. 
The respective softening lengths are 0.3 and 1.2 h~ kpc, and 
the minimum galaxy baryonic mass resolutions are 2.1 x 10 6 
and 1.4 x 10 8 Mq. We will show that the primary enrichers 
of the IGM are galaxies with masses ~ 10 7_8 ' 5 M Q , which 
are statistically sampled best in the 16/i _1 Mpc box. 

The dl6n512vzw simulation was run to z = 3, however 
we only consider outputs down to z — 5 here. This simu- 
lation was run on the Intel 64-bit Abe Cluster at the Na- 
tional Center for Supercomputing Applications on 128 pro- 
cessors, requiring 32,500 processor hours. The d8n512vzw 
and d32n512vzw simulations were run only to z = 5.5 
and 4.5 respectively, each consuming approximately half the 
above processor hours. We run the Spline Kernel Interpola- 
tive DENMAX (SKIDfl group finder on the z =8, 7, 6, and 
5 outputs to identify galaxies. 



3 PHYSICAL PROPERTIES 

In this section we explore key physical and observational 
predictions of our new simulations and relate them to find- 

3 http:/ /www-hpec. astro. washington.edu/tools/skid. html 



ings from our previous lower-resolution simulations. This is 
intended to demonstrate continuity with our previous works, 
as well as to provide a more global context for our explo- 
ration of high-z metal-line absorbers. We begin by present- 
ing the z — 6 galaxy stellar mass (M*) function in our 
two simulation boxes, and relating this to earlier papers 
where we fit observables best with a momentum-driven wind 
model. We show that the 16 /i -1 Mpc box is the best one for 
studying high-z IGM metals, despite the fact that it is not 
large enough to sample the space density of the most lumi- 
nous z ~ 6 galaxies. Finally, we show global evolutionary 
trends in star formation, wind properties, and metal enrich- 
ment. 

3.1 Galaxy Mass Function 

In DFO06 we explored the properties of reionization-epoch 
galaxies (« = 9 — > 6), finding that the mzw wind 
model most capable o f reproducing observations (specifically 
iBouwens et al.l 2006) compared to the constant wind (cw) 
and no wind (nw) cases. The mzw model is similar to the 
vzw model, the only difference being that the wind speed at 
a given a is constant in mzw while it has a random spread in 
vzw. The key point is that r\ varies with a the same way in 
both models, and as we showed in DFO06, r\ governs high-z 
SF by removing gas that would othe rwise form into stars. 
The co nstant wind model is that of ISpringel fe Hernquistl 
|2003bl ), with w wind = 484 km/s and 77 = 2. 

We plot the stellar mass function of z = 6 galax- 
ies ($(A/»)) in Figure Q] for the d8n512vzw (orange), 
dl6n512vzw (black), and d32n512vzw (green) simulations. 
We compare it to the w8n256vzw simulation (dot-dashed 
magenta), which is part of the w-series simulations ex- 
plored in DFO06. The w8n256vzw simulation has nearly 
the same resolution as the d!6n 512vzw simulat i on, b ut uses 
the first-year WMAP results dTegmark et al . 2004). This 
cosmology (fi = 0.3, fi A = 0.7, fi b = 0.04, H = 70 
km s _1 Mpc~\ (78 = 0.90) produces almost double the 
stars over M* = io 70-8,5 Mq owing to its higher ag. In 
the dl6n512vzw simulation, we find 1.83% of baryons in 
galaxies and 0.13% in stars; the respective numbers in the 
w8n256vzw simulation are 3.10% and 0.26%. 

We use $(M«) to argue that we are resolving the galax- 
ies responsible for metal enrichment. Stellar mass can be 
used as a direct proxy for the integrated metal production if 
two conditions are met: (i) All metals (those not trapped in 
stellar remnants) are produced in short-lived, massive stars 



Tracing the Reionization Era with Metal Lines 5 



9 



10 



T 



I i i i | i i i 
dBn512vzw 
dl6n512vzw 
d32n512vzw 
w8n256vzw 




7 8 9 

log[M. (M )] 



10 



Figure 1. The 2 = 6 stellar mass function of galaxies in 
the d8n512vzw (orange line), dl6n512vzw (black line), and 
d32n512vzw (green line) boxes in units of #/log(c(Af»)/Mpc 3 
(comoving). Corresponding dotted vertical lines show the galaxy 
mass resolution limit. For comparison we show the mass function 
of our older w8n256vzw model (dot-dashed magenta line, from 
DFO06), which is higher due to the different cosmology. The solid 
red line is a least-squares fit to our three simulations (Equation 
0}. We define the characteristic M, as where the fit matches the 
slope of a $ oc M^ 1 power law, which is at M* = 10 7 ' 7 Mq. The 
bottom panel shows a histogram of the stellar mass. 72% of the 
mass is in galaxies with M* = 10 6,7— 8 ' 7 Mq. 



(i.e. Type II SNe), which is a fair assumption given that 
longer lived stars (AGB stars) have not had time to enrich 
much by z ~ 5; and (ii) stellar mass is a direct proxy for 
the integrated number of Type II SNe. The stellar mass con- 
tained in Gadget-2 star particles consists only of long-lived, 
non-SN stars, so the mass in star particles is a good record 
of the integrated Type II SNe contribution. 

We find that a two-dimensional polynomial fit to the 
logarithmic stellar mass function provides a good, simple 
fit: 

log[$(M„)] = -8.70 + 2.951ogM„ - 0.255(logM„) 2 . (3) 

Shown as the red line in Figure [T] this fit obtains a slope 
of —1 at M* = 10 7 ' 7 Mq , making this the characteristic 
M, at z ~ 6 where the most mass per logarithmic M* bin 
is held. We cannot fit the mass function over 4 dex effec- 
tively with a Schechter function, because the high mass end 
does not display an exponential cutoff and the low-mass end 
shows a gradual t urnover owing to filtering by the cosm ic UV 
background (e.g. lThoul fc Weinberdll99a : lllievll2007h . Inte- 
grating Equation |31 we find that 72% of mass is in galaxies 
with M* = io 6 ' 7 " 8 ' 7 M . To illustrate the stellar mass dis- 
tribution, the bottom panel of Figure [T] shows the fractional 



mass in each bin. The median M* in terms of integrated 
mass is M* = 10 7 ' 6 Mq for the dl6n512vzw simulation, and 
M* = 10 7 ' 7 M when including all three simulations. The 
d32n512vzw box is better able to account for the cosmic vari- 
ance of more massive galaxies, but its mass resolution is com- 
parable to the median stellar mass, so it cannot account for 
early enrichment from small galaxies. The d8n512vzw run, 
conversely, does not produce many galaxies above 10 8 Mq, 
so it does not sample the high-mass end sufficiently. The 
dl6n512vzw run is well-suited to resolve the relevant galaxy 
mass range, so we employ this run from here on in. 

3.2 Star Formation Rate Function 

A high-z galaxy's SFR is more easily measured than its stel- 
lar mass, since the former can be estimated from optical data 
while the latter requires near-infrared data. Our simulations 
predict a tight relation between SFR and M* at z = 6, where 



SFR = 9.2 M» M Gyr" 



(4) 



This nearly linear relation is similarly tight in our sim- 
ulations at all redshifts (DFO 06, see their Figure 4; 
iFinlator et~al] l200fj ; iDavel I2008T ), with a low scatte r, and 
is mostly insensitive to outflow model (|Dave1 [2008). This 
relation arises from the dominance of smooth, filamentary 
cold mode gas accretion, where the SFR is regulated by the 
dynamical ti mescale for gas in fall in slowly-growing halo po- 
ten tials (e.g.lKeres et al ■1l2005l ); therefore the SF law we scale 
to (|Kennicutt 19981 ) has little effect on the SFR, chemical 
production, and feedback. Mergers are a sub-dominant pro- 
cess for fueling galaxies, particula rly at early times and in 
smaller galaxies (|Keres et al.| [2008). which is why the scatter 
is small. 

In Figure [5] we show our dl6n512vzw simulated rest- 
frame 1350 A UV luminosity function at z = 6, which is 
effectively a SFR function (including e xtinc tion). We com- 
pare t o data from lBouwens et al.l l|2006l ) and lBouwens et all 
|2007| ). We compute simulated g alaxy luminosities by pass - 
ing their SF histories through the lBruzual fc Chariot! (|2003h 
population synthesis model, and applying a tophat filter 
around 1350A. Dust extinction is critical in this compari- 
son; we apply two extinction laws to our simulated galax- 
ies: (i) the solid line represents a constant dust extinction, 
where E(B - V) = 0.069 , which leads t o 50% flux attenua- 
tion at 1600 A using the ICalzetti et al.l (|2000l ) law (we will 
use this value for our local ionizing field); and (ii) the dashed 
line uses the metallicity-depende nt extinction law calib rated 
from local galaxies introduced in IFinlator et al] l|2006l ) and 
employed in DFO06, which predicts less attenuation for the 
lower mass galaxies. In both cases the agreement with data 
is reasonable. This shows that the star formation rates, and 
hence metal production rates, in our simulated galaxies are 
plausible. 

Numbers across the bottom of this figure show the 
median log[M«] corresponding to a given rest-UV magni- 
tude. The characteristic M» from Equation [3] corresponds 
to ^ 1350 ^ = —17.1, which is roughly one magnitude be- 
low current detection limits. The green dotted lines indicate 
the magnitude range of the parent galaxies of z = 6 C IV 
absorbers, which we identify in H6.3I Hence the originating 
galaxies of currently observed high- z C IV systems are mostly 
below today's detection thresholds. 
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Figure 2. The simulated 2 = 6 rest-UV luminosity function 
(1350A magnitudes AB) from our dl6n512vzw run. Solid line as- 
sumes u niform dust ext i nction of E(B — V) = 0.069, dashed line 
uses the iFinlator et al. I d2006h metallicity-dep endent extinction . 
These are compared to observations from | Bouwens et al.l (120061 *1 
(blue triangles) and iBouwens et al. I ||2007F (magenta squares). 
The agreement is good in the overlap region. The majority of 
SF remains below current detection limits, which implies that 
the sources of high-2 IGM metals have yet to be detected. The 
range of M f° r parent galaxies of 2 = 6 C IV absorbers are 

demarcated by the dotted green lines. Median logarithmic stellar 
masses at a given magnitude are listed along the bottom. 



3.3 Cosmic Star Formation 

Figure |3 (top panel) shows evolution of the cosmic SFR den- 
sity (SFRD) per comoving volume, indicating a rise of a fac- 
tor of three over the range z = 8 — > 5. The values are in the 
range of data, though the observations span a large range, in 
part owing to differences i n the adopted limiting luminosity 
in each case. For example. IBouwens et al.1 l|2007ll finds 6.7x 
more SFRD when they reduce the limiting luminosity by a 
factor 7.5x for the z-band dropouts ((2) = 7.4). Most of the 
z — 6 stellar mass and SF in our simulations remain below 
the detection limits of current surveys. At z = 6, we find 31% 
of the SFRD arises from galaxies with M« > 10 s M© , 39% 
from M* = 10 7-8 Mg, and the remaining 30% from masses 
below this. At z = 8, 61% of the S FRD arises in galaxies 
with M» < 10 7 M Q . We note that IBouwens et al] (120071 ) 
examined the evolution of the cosmic SFRD and found that 
our momentum-driven wind model provided the best match 
to the observed evolution from z ~ 7 — * 4, of the models 
considered there. 
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Figure 3. Top: The cosmic SFRD from all galaxies in the 
dl6n512vzw simulation (solid line), compared to data assum- 
ing a Chabri er IMF at high r edshift. Data are f rom the 
i HopkinsI [120041) com pilation (X's ). iFontana et al.l (|2003|) (circle). 



I Bunker et al 



(120041) (triangle), Thompson et al.l d2006h (pen- 
tagons). and lBouwens et alj feoOS ) (squares). Average wind prop- 
erties (solid lines) in the next three panels show typical t> w i n d = 
200-250 kms" 1 with mass loading factors r\ ~ 6 — 7. The energy 
in outflows is rj 20% of the available supernova energy at these 
redshifts, assuming a Chabrier IMF. One sigma dispersions are 
indicated by dashed lines. The baryonic mass of galaxies respon- 
sible for winds (bottom panel) increase 6x between 2 = 10 — > 4.5 
(bottom panel). 



3.4 Winds 

The adopted dependence of momentum-driven wind scalings 
from galaxy mass results in the median wind speed growing 
and median mass loading factor declining towards lower red- 
shifts (OD08). This has direct ramifications for reionization- 
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era galaxies, since SFRs are attenuated by a factor of (1 + 7?) 
relative to the no- wind case. This is because -r-7— of infalling 
gas that would otherwise be converted into stars, instead 
is ejected in outflows. This factor assumes that outflows do 
not return to galaxies, which is a fair assumption at high 
redshift because outflow material generally does not recy- 
cle back into galaxies within a Hubble time at these epochs 
(OD08). 

The lower panels of Figure [3] show that high-z wind 
properties have a relatively weak redshift dependence from 
z = 10 — » 5: Wwiad increases by 20% (second panel) and 
r\ declines by 17% (third panel), despite an increase of 6x 
in the average baryonic mass of a galaxy driving the winds 
(bottom panel, i.e. weighted by M m i n d). The reason is that 
a is proportional to M ga i x H(z), and the declining H(z) 
turns out to roughly cancel the increase in {M ga i) during 
this epoch, resulting in only modest evolution of wind prop- 
erties. Fiducial values during reionization in our model are 
therefore u w ind ~ 200 kins -1 and r\ ~ 6. The energy bud- 
get of momentum driven winds remains roughly 20% of the 
total SN energy assuming a Chabrier IMF (fourth panel). 
High-z winds require substantial energy input, but the wind 
speeds are relatively mil d compared to low- z outflows from 
LIRGs and ULIRGs (e.g. lMartinll2005al ; lRupke et al.ll2005h . 
As before, we find that substantial mass ejection is a key 
feature of our wind model: the amount of material ejected 
at these epochs is many times that which forms into stars. 



3.5 Metallicity Distribution 

One reason why enriching the IGM is easier at high-z is that 
the physical scale of the Universe is small. OD08 showed that 
in our wind model, winds propagate a similar physical dis- 
tance (~ 60 — 100 kpc) for all redshifts and galaxy masses; 
the comoving volume a galaxy can enrich therefore goes as 
(1 + z) 3 . Figure [3] illustrates the rapid growth of metallicity 
at IGM overdensities from z = 8 — ► 5. The gradually ris- 
ing Z — p trend in place at all four redshifts is one of the 
defining signatures of the momentum-driven wind enrich- 
ment. In contrast, a wind m odel with constant u w ind (e.g. 
ISpringel fc Hernquistl2003bl ) will create a flatter metallicity- 
density relationship with either a sharp drop-off at lower 
overdensity (cf. Figure 10 of OD06) or a rise in metallicity 
if sup erwinds are even stronger (cf. Figure 8 of ICen et al.l 
(|2005|)) signifying efficient enrichment of voids. Under our 
model the IGM metallicity grows between z = 8 — > 5 by a 
factor of 7x at an overdensity (5 = p/p — 1) of 80 and 15x 
at 8 = 8. 

Despite some metals reaching cosmic mean densities, 
most of the volume remains unenriched, especially the voids. 
The cumulative volume filling factors of metals in the bot- 
tom panel of Figure [4] show this, as even regions with 
Z > 10 -3 Zq fill only about 1% of the simulation vol- 
ume at z = 5. Filling factors increase only modestly to 
lower z (OD08), which is somewhat remarkable given that 
these models can match most available IGM metal-line data. 
Hence the presence of diffuse IGM metal absorbers does not 
necessarily indicate widespread metal distribution. We will 
return to this point, and the inhomogeneous enrichment it 
implies, later. 

Regions of galactic overdensities (p/p > 10 3 ) are gener- 
ally enriched rapidly to Z ~ 0.1 Zq or higher. The metal- 
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Figure 4. Top: The average metallicity-density relationship at 
z = 8, 7, 6, 5 from the dl6n512vzw simulation. Z is the total mass 
fraction of metals; Zq = 0.0122. A clear gradient is established 
early on, and as time passes the diffuse IGM is enriched more 
quickly than high-density regions. Bottom: The cumulative vol- 
ume filling factor of metals versus metallicity, indicating an inho- 
mogeneous distribution of metals with only 1% of the Universe 
enriched to 10 — 3 Zq by z = 6, despite some metals spreading to 
near the cosmic mean density. 



licity does not grow much at such overdensities, which is 
reflected in the lack of predicted evolution in the mass- 
metallicity relation from z = 8 — + 6 in DFO08 (see their Fig- 
ure 5). The fraction of cosmic volume enriched to Z — Zq 
is 0.01-0.02%, which is actually a factor of 1.5-3 x higher 
than that at z = 0.25 (|Oppenheimer fc Pavel lioogh . It is 
astounding to consider that the volume fraction of the Uni- 
verse enriched to around solar metallicity may change very 
little between z = 8 — » 0; the key difference is the cor- 
responding regions at low redshift contain orders of mag- 
nitude more baryons. This is a result of our moment um- 
driven wind model based on locally observed outflows. 
iDave fc Oppenheimerl l|2007h showed these winds efficiently 
enrich IGM overdensities at high-z, some to ^ Zq, while by 
z — most of these metals recycle back into galactic envi- 
ronments and later winds often cannot escape out of galactic 
haloes (OD08). 

Many of the trends noted above can be seen in Figure [5] 
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which shows the locations of galaxies in the upper left and 
the metallicity-density distribution in the upper right in a 
25 kms~ 1 -thick slice extending across the dl6n512vzw sim- 
ulation box. Gray indicates zero metallicity in these frames, 
which corresponds to the vast majority of the IGM, and 
shows the nascent cosmic web. Most of the volume remains 
unenriched at these redshifts, while regions within early 
galaxies become enriched up to a tenth solar or more. Clus- 
tering between early galaxies and enriched portions of the 
IGM is strong. This illustrates that the galaxy-absorber con- 
nection at high-z is extremely tight, as we investigate further 
in SJO] 



4 IONIZATION BACKGROUNDS 

To make predictions for observable metal-line absorbers, 
we must obtain ionization corrections for relevant metal 
ions. This requires that we have a model for the photo- 
ionization rate at each point in the simulation volume. Ide- 
ally, we would like to self-consistently track the local ioniza- 
tion spectrum in these simulations. But this must await a 
full time- dependent, non-equilibrium radiative transfer im- 
plementation, which we a re in the process of developing 
i|Finlator et all 120081 . [20091 ) but do not yet have. Hence for 
this work, we will consider three ad hoc cases for the ioniza- 
tion background (dubbed No Field, HM2001, and Bubble) 
that span a plausible range. 

Not e that our simulations are originally run with a 
uniform lHaardt fc Madaul |200ll ) ionizing background, and 
we are only imposing these backgrounds in post-processing. 
This obviously does not properly capture the dynamical ef- 
fects of different pressure forces in our three cases. How- 
ever, since dark matter is dynamically dominant in the rela- 
tively low-density regions giving rise to IGM absorbers, post- 
processing the background to obtain ionization corrections 
is a good approximation. 

4.1 Ionization Cases 

No Field: Our first case assumes a fully neutral IGM, 
where no ionizing background exists at wavelengths below 
the Lyman limit, while lower energy photons stream freely. 
This extreme case may be applicable if metal injection 
precedes ionization fronts, or more realistically if the 
IGM recombines in regions enriched and ionized at earlier 
times. Oxygen remains in its ground state due to the 
Ol ionization potential (13.618 eV) being barely higher 
than the Hi ionization potential, while Cl (11.260 eV) 
and Si I (8.151 eV) are completely ionized to the second level. 

HM2001: Our second case assumes reionization is complete 
and a spat ially-uniform bac k groun d applies everywhere. 
We use the lHaardt fc Madaul (|200ll ) background assuming 
a contribution from quasars and 10% of photons above 
the Lyman limit from star forming galaxies generated by 
their Cosmic Ultraviolet Background (CUBA) program; at 
z > 6, star-forming galaxies dominate the ionizing flux. We 
divide the total flux by a factor of 1.6 as we did in OD06 
in order to match the summed flux decrement of the Lya 
forest down to z = 2. An assumption of a spatially-uniform 
ionization background is probably far from correct at 



z > 6, but this is intended to represent another extreme 
case like No Field. Ionization fraction s for metal lines are 
calculated using CLOUDY-generated (|Ferland et al.lll998t ) 
ionization tables as a function of density and temperature, 
as described in OD06. 

Bubble: The last case assumes that the ionizing flux of the 
nearest galaxy dominates, which may be the most realistic 
scenario during reionization. This is applicable if metals re- 
main near enough to their parent galaxies to be dominated 
by the local ionizing intensity, which we will show is often 
true since metals typically reside within 30 kpc of the orig- 
inating galaxy. Of course, detailed radiative transfer effects 
of escaping ionization and subsequent recombinations will 
play a critical role, so this scenario should be considered a 
crude proxy for a patchy reionization model. 

To implement the Bubble field, we consider the ionizing 
field from only the nearest galaxy, weighted by the inverse 
square law. We as sume an average SEP per stellar mass by 
taking the average [Bruzual fc Chariot] (|2003l ) stellar synthe- 
sis spectrum of all the stars in the simulation box for a given 
redshift, accounting for their ages and metallicities. We then 
calculate the impinging flux using the stellar mass and dis- 
tance of the nearest galaxy, which is identified for each SPH 
particle. Note that the same ionizing flux per unit stellar 
mass is assumed to emanate from each galaxy at a given 
epoch; we do not account for variations in nearest galaxy 
SEDs, which is another approximation required to make the 
problem computationally tractable. 

We attenuate the emission from galaxies by dust extinc- 
tion. The assumed extinction has a fixed value of E(B—V) = 
0.069, as shown in Figure [2] to broadly match the observed 
rest-UV LF. To calculate how the emergent SED is atten- 
uated, we use the theoretic al dust absorption cur v e from 
the X-ray to the far-IR of IWeingartner fc Dralnel (|200lh . 
using their grain size distribution of silicate and carbona- 
ceous grains found to fit the Small Magellanic Cloud (SMC) 
bar. We use this curve because z ~ 3 LBGs exhibit the 
smaller grain size distributions of SMC-like dust (Viih et al.l 
2003) . Note that w e cannot directly apply the empirical 
Calzetti et all (2000) extinction law because it extends only 
down to 1200 A, while we are especially interested in the 
attenuation at energies where Civ and Si IV are ionized (i.e. 
< 400 A). For each SPH particle, the impinging flux is then 
the dust-attenuated emission divided by the distance from 
the nearest galaxy squared. We then use CLOUDY to gen- 
erate the resulting metal ionization fractions, which is then 
used to obtain the optical depth for any ion along the LOS 
as described in OD06. 

We plot the HM2001 field at z = 8 and 6 in Figure El 
along with the Bubble field, where we assume a M» = 1O 7 M0 
galaxy at 10 kpc (typical values) for concreteness. The dif- 
ference between the solid and dotted orange lines for the 
Bubble case at z — 6 illustrates the effects of dust attenu- 
ation. Dust absorption in all of the [Weingartner fc Drainei 
l|200ll ) models sharply peaks at 600-900A, causing a dip. 
There is a clear Lyman break, which in our cases arises 
mostly from stellar atmospheres, and secondarily from dust. 
We note t hat this break is not nearly as strong as that ob- 
served bv lShaplev et all l|2006h in Lyman break galaxies at 
z ~ 3, perhaps indicating that there is a fairly small escape 
fraction of ionizing photons; for instance, lHaardt fc Madaul 
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Figure 5. A 25 krns" 1 slice of IGM spanning the 16 comoving h~ 1 Mpcs of the dl6n512vzw simulation at z = 6. Upper left: The 
locations of galaxies are shown as colored points, where the color corresponds to stellar mass, overlaid on the gas overdensity represented 
in greyscale. Upper right: Color indicates the enrichment level; the vast maj ority of the IG M remains unenriched, and metals are close 
to galaxies. The scale bar indicates the arcsecond scale on the sky using the IWrierrd d2006l) on- line c alculator. Bottom panels: C II (left) 
and Civ (right) column densities in this 25 km/s slice, assuming a uniform lHaardt fc Madaul ll200lh ionization background. Cn traces 
denser regions compared to C IV, which traces the volume distribution of metals extremely well. 



l|200ll ) assume a global 10% escape fraction. However, be- 
cause these data are from lower- z and higher- M» galaxies, it 
may not be applicable here. For simplicity, we do not assume 
an overall escape fraction. 

Comparing the Bubble and HM2001 fields at z = 6, the 
strength of the Lyman break is much less in the Bubble case. 
This is mostly because our typical gala xy is younger and 
more metal-poor than that assumed by lHaardt fc Madaul 
(|200ll ) (0.2 Zq, 0.5 Gyrs), which leads to a harder ionizing 
source. Two secondary reasons for greater hardness are: (i) 
the Bubble Lyman-break attenuation due to dust is less than 
the HM2001 attenuation applied (~ 3 vs. an overall 10% 
escape fraction); and (ii) the dust attenuation declines at 
the ionization potentials of C IV and Si IV (shown) . This leads 
to greater ionization fractions at these higher states in the 
Bubble field, which has observable implications as we will 
show throughout Sj5] 



4.2 Ionization Fraction Behavior 

Before delving into the observables, we present here a primer 
for understanding metal-line observations resulting from our 
two assumed ionization backgrounds, the HM2001 and the 
Bubble fields. The No Field case is straightforward (i.e. ion- 
ization fractions of unity for C II, O I, and Si n), but the often 
nuanced behavior of the other two fields has some perhaps 
unexpected consequences. We consider metals in p — T phase 
space exploring both low and high-ionization species at two 
redshifts that span the redshift range we consider. 

Figure [7] plots colored contours corresponding to the 
HM2001 (left panels) and the Bubble (right panels) back- 
grounds in p — T phase space for C II (solid blue) and C IV 
(dot-dashed red) at z = 8 and z = 5 (upper and lower pan- 
els, respectively). Each contour step represents an increase 
of 20% in the ionization fraction for the corresponding ion. 
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Figure 6. Ionizing background spectra: The HM2001 (dashed) 
ionization background at z = 8 (magenta) and z = 6 (orange) 
are compared to the Bubble field arising from a z = 6 Af* = 
10 7 M@ galax y at 10 kpc. The thin solid orange line represents the 
unattenuated [Bruzual fc Charlotl j2003h -computed SED at z = 8 
to compare with the dust-attenuated (thick dotted orange line) 
SED. For reference, ionization potentials for carbon, oxygen, and 
silicon are indicated, and the relative cross-sections from Hi and 
Hen as a function of wavelength are shown along the bottom. 



Overlaid are thick gray logarithmic metal density contours 
(darker contours indicate higher metal densities) at the re- 
spective redshifts. Overlapping color and gray contours rep- 
resent regions of phase space where absorption in that ion 
will arise for a given ionization field. 

The HM2001 contours are the same at all spatial lo- 
cations within the simulation box at a given redshift, so 
the resulting behavior is simpler to understand. The cor- 
responding contours for both C II and C IV move to higher 
overdensity by almost a factor of 10 from z — 8 — > 5, ow- 
ing to the ionization parameter U increasing. U is defined 
in this paper as the number density of photons capable of 
ionizing to a particular state (e.g. ionizing C ill to C iv) di- 
vided by the number density of a particular atomic species 
(carbon in this case) . U increases ~ x 10 because the phys- 
ical densities decline by a factor of 3.4 from Hubble expan- 
sion, and the ionization intensity increases by ~ x3 be- 
tween z = 8 — > 5. C II evolves to trace a smaller fraction 
of phase space at higher densities; however the metal den- 
sity is greater at these higher densities and at later times, so 
the total amount of C II is expected to be fairly stable from 
z = 8 — ■> 5. The contours extending to higher densities at 
T ~ 10 5 K in Figure [7] correspond to Civ in collisional ion- 
ized equilibrium (CIE). This overlaps few metals, because 
C iv is an efficient coolant and SPH particles evolve rapidly 
through this region, much li ke we find for O VI at z < 0.5 
|Oppenheimer fc Dave|[2009l ). 

As most IGM metals remain near T — 10 4 K, we plot 
the Cn/C and Civ/C ionization fractions (/(Cn), /(Civ)) 
at this temperature in Figure [8] (top panels) , showing the 
HM2001 curves moving upward in overdensity by the factor 



of 10 from z — 8 — ► 5. C iv should dramatically increase as 
the 2 = 8 field is too weak to ionize C iv, while by z — 5 
the gray and C iv contours move to overlap as more met- 
als enrich lower overdensities while U increases. A second 
effect illustrated in Figure [8] is that the maximum /(Civ) 
grows as a result of the harder background with more quasar 
contribution at z — 5 relative to z = 8; more photons ex- 
ist at the ionization potential of C ill (47.9 eV or 259 A) 
to create Civ relative to the Ci potential (11.3 eV or 1101 
A) to make C II, resulting in a larger density range where 
C iv resides. We define and plot in Figure [8] (bottom panel) 
the ratio J ° I J ° , which we call the "carbon spectral 

259A' 1101A 

ratio" . This has significant implications for C n and C iv ob- 
servations. The carbon spectral ratio is 200 x greater for the 
Bubble versus the HM2001 field at z = 8. 

We now focus on the Bubble field behavior. The example 
fields signified by the right panels in Figure [7] demonstrate 
the expected ionization in p — T phase space at 10 physi- 
cal kpc from a 10 7 M© galaxy. These contours shift right 
as M* increases and left as d ga i increases due to the inverse 
square law; e.g. the contours are equivalent to the ioniza- 
tion field from a 10 9 Mq galaxy at 100 kpc. We choose this 
field strength because this is the representative of where our 
simulated metals typically reside. Like the HM2001 field, 
our implementation of the Bubble field depends on only two 
variables: U and T. The difference is that U depends on a 
combination of density, stellar mass, and galaxy distance. 

If we consider the Bubble case shown in Figure [7] at 
2 = 8, there exists moderate overlap with the metal con- 
tours between 1 — 10 x mean overdensity, which are reason- 
able densities at 10 physical kpc from our simulated 10 7 Mq 
galaxies. But note that the calculated incidence of intersect- 
ing at least this close to one of 332 M* ^ 10 7 Mq galaxies in 
the dl6n256vzw box is 0.22A« _1 , so a high-z quasar spec- 
trum would be lucky to intersect one such region. The con- 
tours move to higher overdensity for larger galaxies (or at 
closer distances), but these are even less likely. Meanwhile, 
by z — 5, the space density of M» ^ 10 7 Mq galaxies is x8 
higher and more metals have enriched IGM overdensities 
while the ionization contours have not evolved too much. 
This is why we obtain an increase in C iv absorption from 
z = 8 — ^ 5, as we show in >15.2.2I 

The greater carbon spectral ratio in the Bubble field 
leads to greater /(C iv) at photo-ionized temperatures, as 
shown in Figure [5] At T = 10 4 K, /(C iv) maximizes at 
0.67 at z = 8 and extends over 3.2 decades of density where 
/(Civ) > 0.1, compared to /(Civ) = 0.27 and 1.6 decades 
with /(Civ) > 0.1 for the HM2001 field. The Bubble hard- 
ness declines as stars age and grow more metal rich by z — 5, 
(maximum /(Civ) = 0.57), and rises for HM2001 as the 
quasar density increases (/(C iv) = 0.37, cf. the carbon spec- 
tral ratio evolution in the bottom of Figure [Sj . The greater 
Bubble hardness allows significant ly more Civ, e specially at 
the highest redshifts we explore. ISimcoel l|2006h also finds 
/(Civ) increases at z > 4.5 at p/p = 10 when the dust at- 
tenuated flux from a 10 s Mq galaxy at 100 kpc is included. 
In the end, this results in less evolution for C iv absorp- 
tion in the Bubble model relative to HM2001, as we show 
in £15.2.21 Our Bubble evolution is greater than the ISimcoel 
(2006) example, because the metals are 5-50 kpc from ion- 
izing galaxies. 

A key point applicable to all metal-line IGM observa- 
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Figure 7. Metal phase space plots: Colored contours (solid blue for C II and red dot-dashed for C iv) correspond to the ionization 
fractions, while thick gray contours indi cate the metal dist r ibutio n in phase space. Overlapping regions give rise to absorption in that 
ion. The left panels assume the uniform lHaardt &: Madaul fcOOlh background, and the right panels use the local ionizing Bubble field 
arising 10 kpc from a 10 7 Mq galaxy; the contour spacings correspond to ionization fraction steps of 20% (i.e. /(Civ) =0.2, 0.4, etc.). 
The gray metallicity contours are in logarithmic steps of 0.5 dex. At z = 8 (top panels), Cn should trace metals better than Civ for the 
HM2001 field, while by z = 5 (bottom panels) C IV should become a better tracer of metals in the diffuse IGM. The Bubble contours 
depends on galaxy size and distance: Contours move right as Af* increases and left as d ga i increases (i.e. U tx M*d~^). Civ is always a 
good tracer of IGM metals from z = 8 — > 5 in the Bubble case. 



tions is that while it is well known that taking ratios of 
different photo-ionized specie s can determine spec t ral ha rd- 
ness (e.g. ISchave et ail 120031 ; lAgiiirre et all |2004 120081 ). a 
single species' ionization correction alone can vary dramat- 
ically due to the spectral shape. It is risky therefore to as- 
sign a fiducial ionization correction, especially considering 
we know so little about the ionization background for high- 
ionization species; as a result, an observable such as f2(Clv) 



could easily change by a factor of a few, as we stress in 

mm 

To visualize the spatial distribution of metal ab- 
sorbers, we illustrate the absorption column densities in two- 
dimensional maps of the 25 km s" 1 slices in the bottom pan- 
els of Figures [S] assuming a uniform HM2001 background. 
Cll is confined to overdense regions close to galaxies, while 
C IV is a more capable tracer of the IGM volume enriched, 



12 B. D. Oppenheimer, R. Dave, & K. Finlator 



l 

0.8 

^ 0.4 
0.2 

0.8 
0.6 



-1 



> 

^ 0.4 K 
0.2 




.y 



.7 



■ - I — I — I H - ■■{■■ F - H - 



I I I I 



- z = 8 Local Ionizing Bubble 

z=5 (10 7 M at 10 kpc) 

—.— Haard & Madau 2001 



-./- ' ? - 

■ / 



••Ns 



log [p/p] 




M 
O 



6 7 
Redshift 

Figure 8. Top panels: C II (upper) and C IV (lower) ionization 
fractions as a function of overdensity at T = 10 4 K for the 
HM2001 (dashed lines) and Bubble (dotted lines) at z = 8 (ma- 
genta) and z = 5 (green). The HM2001 field evolves more as 
quasars provide increasing intensity above 47.9 eV at later times. 
The Bubble field shows less evolution, and generally C IV traces 
lower overdensities. Lower panel: The mean carbon spectral ra- 
tio = J ° I J ° for our two ionization fields as a function of 
259A' 1101A 

redshift. The star forming galaxy-dominated Bubble field is sig- 
nificantly harder than the HM2001 field with quasars, because 
HM2001 is dominated by galaxies at high- z and the attenuation 
above 47.9 eV is greater than the Bubble attenuation due to dust. 
This results in overall lower values for /(C iv) in the HM2001 



albeit mostly at column densities that are below current ob- 
servational limits. The scarcity of pink in the bottom right 
fra me indicating column den sities similar to thos e observed 
by iRvan- Weber et all (|2006l ) and ISimcoei jjool ) indicates 
just how rare detectable C IV is. 

We demonstrate evolution using a zoom-in to the second 
most massive galaxy in our box in Figure [9] This visualiza- 
tion shows the evolution of the surrounding 4 x 4/i~ 1 Mpc x 
25 kms -1 region using the same scales as Figure [S] A clear 
bipolar outflow emanates perpendicularly to a filament feed- 
ing the growth of a M„ = 10 9 Mq galaxy (at z = 6). 



Strong Cn absorption traces the overdense regions around 
the galaxy at all epochs, while C IV evolves to trace the ex- 
tended IGM enrichment; observed C IV evolves faster than 
C II between z = 8 — » 5. The smaller galaxies in the bottom 
right of this frame, M* ~ 10 7-8 Mq, are the most prolific 
IGM enrichers at high-z due to their large numbers and high 
mass loading factors. 



5 OBSERVATIONAL PREDICTIONS 

Observations of the high- 2 IGM are leaping forward with the 
advent of near-IR spectroscopy, along with the discovery of 
quasars at z > 6, with an eye to th e tantalizing possi bil- 
ity of using gamma-ray bursts (GRBs: lKawai et aLll2006l ) to 
briefly illuminate IGM absorbers at very high redshifts. We 
present here predictions for high-z metals line transitions 
occurring redward of Lya, where the quasar continuum is 
free of Hi absorption. We consider five species already ob- 
served at z 5 with rest wavelengths (Ao) redward of Lya: 
Cn (1334.5A), Civ (1548.2, 1550.8 A), Oi (1302.2 A), Sin 
(1260.4 A), and Siiv (1393.8, 1402,8 A) at our four chosen 
redshifts: z = 8, 7, 6, and 5. Three ionization fields intro- 
duced in H4.1l are considered: No Field, HM2001, and Bubble. 

We use our quasar absorption line spectral generator 
specexbin (see §2.5 of OD06 for details) to produce con- 
tinuous angled LOSs through our simulation volume, each 
covering Az — 0.9 for each of the four redshifts. We make 
an "ideal" -quality sample of 70 LOSs, assuming a high S/N 
ratio of 50 for each R — 50,000 resolution element, and 
a moderate instrumental resolution of 15 kms . This is 
comparable to what is achievable on current optical spec- 
trographs on 8 — 10 meter class telescopes, which can trace 
Cn out to z ~ 6.5 (e.g. BSRS). At higher redshifts, near- 
IR data will be required. The current generation of near-IR 
spectrometers have resolutions of R = 5000 — 10, 000 and 
a S/N = 10 — 20 per resolution element for z ~ 6 quasars 
(R. Simcoe, private communication). The next generation 
should achieve R — 20,000 with much higher S/N using 
20-30 meter class telescopes, e.g. the near-Infrared Multi- 
Object Spectr ograph (IRMOS) on th e Thirty Meter Tele- 
scope (TMT) l|Eikenberrv et al.ll2006l ). Our simulated spec- 
tra can be used both to compare to current optical spectra, 
as well as make predictions for future instruments. 

We first consider column density distributions for three 
of the five ionization speci es. We fit Voigt p rofiles using the 
automated fitter AutoVP |Dave et al.lll997l ). generally with 
little uncertainty since metal lines are typically unsaturated 
and unblended. Then we sum column densities into a cosmic 
mass density as a function of redshift for each species. We 
present predictions for equivalent widths next, for compar- 
isons with large samples of z > 6 metal absorbers that may 
obtained in the near future. We then consider some sim- 
ple absorber ratios, emphasizing what can be inferred about 
the ionizing conditions. Lastly, we show some sample ab- 
sorption line profiles and explain the velocity profiles of sys- 
tems. Along the way, we compare to published z > 5 metal 
absorber data, and offer some interpretations based on our 
models. We consider the observed completeness limits as- 
suming S/N=10 and 20 per resolution element at R = 5000, 
and we confirm that the results we show in our high quality 
spectra hold at lower resolution and S/N. 
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Figure 9. The evolution of a 4 X 4ft x Mpc X 25 kms 1 slice, centered on a galaxy that grows from 2.4 X 10 8 M to 1.8 X 10 9 Mq from 
z = 8 — > 5. The global metal volume filling factor increases by ~ X 10 from z = 8 — > 5, but this slice shows less growth owing to early 
enrichment from the massive galaxy at the center. C II traces higher overdensities in the proximity of galaxies, while C IV evolves to trace 
the more diffuse IGM by z = 5, including the bipolar outflow. M* = 10 7— 8 Mq galaxies such as those in the bottom right of the frame 
represent the most prolific enrichers of the z = 6 IGM, since they exist in greater numbers and have higher mass loading factors. Visit 
http://luca.as.arizona.edu/~oppen/IGM/hiz.html for simulation movies. 



5.1 Column Density Distributions 

The column density distribution (CDD), d 2 n/dXdN (ion), 
is the most basic metal absorber counting statistic. Here n 
is the number of lines with column density between iV(ion) 
and iV(ion)+dAT(ion), and dX is absorption pathlength. Fol- 
lowing general practice among observers, we sum together 
all individual Voigt profile components within 100 kms -1 
of each other into systems. This minimizes differences due to 
vagaries of Voigt profile fitting, and allows more straightfor- 
ward comparisons among data sets with varying resolution 
and signal-to-noise. For dX we use 



X{z) 



3Qi 



VS!m(1 + z) 3 + !1a. 



(•») 



dX for our adopted Qm = 0.25 cosmology is almost double 
that in an Einstein-deSitter universe (the Ha term is negli- 



gible at these redshifts). We scale other observations to our 
assumed cosmology where necessary. 

Figure [TOl shows CDDs for oxygen and carbon using dif- 
ferent choices for the ionizing field at z = 8 (left) and z = 6 
(right). The first ion considered is Ol in the top panels, 
which has a unique association with Hi becaus e these two 
speci es remain in charge-exchange equilibrium (|Osterbrockl 
119891 1 due to their neut ral ionizat ion p otentials being only 
0.15% apart in energy. lOhl ((2002) and iFurlanetto fc Loebl 
(2003) hypothesize the existence of an "O I forest" tracing 
the neutral metals that have enriched dense IGM regions via 
early enrichment. Higher ionization states easily recombine 
at the high densities of the early Universe. Dozens, even hun- 
dreds of lines with JV(0 1) > 10 14 cm -2 should be detectable 
if the IGM is neutral and enriched with a large filling factor 
of metals, but BSRS find very few O I lines. No O I systems 
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Figure 10. Predicted differential column density distributions for Oi, C II, and Civ systems at z = 8 (left) and z = 6 (right), for our 
three ionization field cases. 50% detection thresholds are indicated by gray lines for S/N = 10 (dashed) and 20 (dotted) per resolution 
element in an R = 5000 spectrum. The two values on the right side of each panel are log[B] and a from Equation \E\ if both numbers 
are negative, or if there are too few lines for a power law fit the line frequencies per Az in the lower and upper column density bins 
if both numbers are positive. Oi (upper panels) shows frequent strong lines and a shallow decline toward high N(Ol) if all oxygen 
is neutral, while significantly less arises in the ionized cases. Cn, middle panels, is assumed to have no ionization correction in the 
No Field case, showing much the same CDDs as neutral O I. C II is more prevalent in the Bubble than O I due to this state's higher 
ionization potential. C IV, bottom panels, shows more strong absorbers for the Bubble field, while the HM2001 uniform field predicts 
much stronger evolution. The long dashed cyan CDD at z = 6 is the case where a metallicity floor of Z = 10" 3 Zq is added, resulting 
in appreciably more lines below 10 13 ' 5 cm -2 . The black square in the bottom r i ght p ane l correspond to the t otal n umber of components 
(3) with AT(Crv) = lo 14 ' 2 " 15 cm -2 in the combined sample of lBecker et al.l 00091) fclRvan- Weber et al.1 l|2009h over a pathlength of 
AX = 27.4. The black upside-down triangle is the upper limits fro m thelBecker et al.l J2009t) null detection at JV(C iv) = 10 13 ' 4 cm -2 
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ab ove 1 15 c m 2 are found either, despite such predictions 
bvlOhl (|2002l ). 

Our No Field case (upper panels of Figure I10[) shows 
that O I lines should be common at z — 8, despite our metal 
filling factors being under 1%. We predict line frequencies of 
1.9 and 0.6 Az" 1 for 10 14 ^ 7V(Ol) < 10 15 cm" 2 and 10 15 > 
N(0 1) respectively at z = 8. The corresponding frequencies 
rise to 5.4 and 1.9 Az by z — 6. Ol provides the clearest 
contrast of any species between a neutral and ionized IGM; 
the Bubble field yields frequencies of 0.4 and 0.3 Az" 1 above 
N(0 i) = 10 14 cm -2 at z = 8 and 6 respectively. This clear 
contrast may already have been observed by BSRS, as all 
of their O I absorbers detected at z > 6 are along one of 
their three LOSs probing these reds hifts. This could suggest 
a patchy reionization scena rio (e.g. iFurlanetto fc Ohll2005l ; 
lllievl 120061 ; iLidz et all 120071 ) where individual LOSs probe 
large-scale structures at different ionization states. 

We fit a power law of the form 



i'N\ 

d n/dXdN "(ion) = B x 



No 



(6) 



to the 1 distribution by summing systems into two bins 



(JV(Oi) 



10 1 



and N(0 1) 



10 



14.45-15.2 



and 



finding the best-fit line. We list log[£>] and a for each CDD 
in Figure [10] if there are at least 10 systems in each bin; oth- 
erwise we list system frequencies per Az of the lower and 
upper bins (positive numbers). We set No = IO 130 cm -2 
for the normalization. The No Field O I z — 6 CDD has 
a — — 1.1 over these column densities, which translates to 
85% as many systems in the upper O I bin compared to the 
lower bin. A power law where a > —2 indicates more metals 
at higher column densities, meaning this power law must 
turn over so that the summed O I is finite. Power law slopes 
can yield information about the nature of the metallicity 
distribution and ionizing source, and provide an easy way 
to compare datasets to simulations. Not only does the Bub- 
ble field create fewer O I lines at every column density, it 
has a much steeper power law, which we cannot calculate 
accurately as there are too few absorbers in the higher bin. 
Oxygen and hydrogen are easily ionized out of their neutral 
states by any ionizing background above 1 Rydberg. 

C 11 CDDs shown in the middle panels of Figure 
1101 recreate similar behavior to O I for the No Field case, 
but has a much higher frequency for the ionization field cases 
since the C I ionization potential lies below 1 Ryd, unlike O I. 
The z = 6 frequency of the lines with 7V(C 11) > 10 140 cm" 2 
is 4.5, 0.8, and 0.4 Az" 1 for the No Field, HM2001, and 
Bubble fields, respectively. Furthermore, the slopes of the 
CDDs between JV(Cn) = 10 13 2 " 14 7 cm" 2 may help con- 
strain the ionization field, being a = —1.1, —1.5 and —1.9, 
respectively. The Bubble field produces fewer strong C 11 lines 
than HM2001, because the metals at z — 6 still reside close 
enough to the local ionizing source that carbon is ionized 
to higher states by the greater intensity above the Lyman 
limit. 

C IV has a higher frequency of lines with N(C iv) ^ 
10 14 ' 2 cm" 2 (0.9 vs. 0.2 Az' 1 ) and a shallower slope (a = 
-1.4 vs. -2.0 over iV(Clv) = io 130 " 14 - 5 ) for the Bubble 
field compared to the HM2001 field in the bottom pan- 
els of Figure 1101 There is no C IV in the No Field case, 
as shocks from structure formation rarely heat the gas to 
temperatures where C IV is collisionally ionized at these 



early epochs. We plot a black square from t he th ree ob- 
served C IV absorbers of I Ryan- Weber et all (|2009T > where 
iV(Clv) = io 14 - " 15 cm" 2 and include the additional path- 
length from t he SP SS J0002+2550 LOS observed only by 
iBecker et all (120091 ) (total AX = 27.4). We find an ob- 
served frequency of 0.5 ± 0.3AZ" 1 for these absorbers av- 
eraging (z) — 5.75. To compare directly to our models, 
we (i) multiply the frequ ency of our C IV absorbers by the 
I Ryan- Weber et aU (|2009T ) sensitivity function (their Figure 
6), and (ii) estimate the z = 5.75 Civ frequency by tak- 
ing a log-weighted average of our z = 6 and 5 model fre- 
quencies, because strong C iv absorbers rapidly increase in 
number from z — 6 — > 5 (especially for HM2001). Our 
JV(CIV) = 10 14 " 15 cm" 2 frequencies are 1.2 and 0.4 
Az" 1 for the Bubble and HM2001 models respectively. The 
uniform field is consistent with the observed frequency, while 
the local ionizing bubble actually over-predicts strong C iv 
absorbers. We actually prefer the the Bubble model, because 
we will argue IGM metals are not completely ionized and 
that a significant fraction of them are in lower ionization 
species, either recombined or shielded from UV flux. 

We also plot the maximum allowed frequency of 
jV(Cl v) = 10 13 ' 4 cm" 2 systems from the null IBecker et all 
(2009) detection (upside-down triangle); this is the column 
density where they calculate they have 50% completeness. 
We appear to over-predict the amount of these weaker ab- 
sorbers in either model, although by how much will have 
to wait for future observations. The long dashed cyan line 
is the case where a 10~ 3 Zq metallicity floor is added to 
the simulation; we discuss this more in H5.3I The fact that 
IBecker et al. I (|2009h finds no weak C iv lines in the largest 
pathlength yet sampled at z > 5.5 indicates a large volume 
of the IGM is not enriched much above Z > 10~ 3 Zp . 

For comparison with lower z, we plot Son gailal (l200ll ) 
power law fit to Civ absorbers at z = 2.90 — 3.54. There 
exists less C iv at all column densities in our simulations, 
which agrees with our reproduction of a declining Q(C iv) 
at z > 5 as we discuss next. 



5.2 Metal Mass and Ion Densities 

We now consider the summed mass densities as a way to 
illustrate the evolutionary differences between ionization 
fields. The evolution of C iv for the two fields differs since the 
weak HM2001 field at z = 8 is incapable of ionizing overden- 
sities containing carbon up to C iv, creating very few strong 
absorbers, while the Bubble field has a harder spectrum early 
on (cf. red contours in upper panels of Figure [7)l . 

An ion's global mass density may be obtained by inte- 
grating the total column density of their systems and divid- 
ing by pathlength, as follows: 



il(ion, z) 



H m spccics £iV(ion, z) 
AX(z) ' 



(7) 



Cpcrit 

where Ho is the Hubble constant, c is the speed of light, 
pcrit is the critical density of the Universe, and m spcc i cs is 
the atomic weight of the given species. 

For comparison to the data, we count systems over 
two decades of column density chosen to correspond to 
the species' detectability. We count lines for C iv over 



iV(Civ) = 10 1 



cm , and scale other species' col- 



umn density range using the ratio of the oscillator strength 
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of a line (g, the strongest of a doublet) times its rest wave- 
length (g x Ao) for the considered species versus that of 
C IV. For example, the range considered for O I works out 
to be JV(Oi) = io 13 - 7-15 ' 7 cm -2 , because the resulting 
Civ : Ol ratio is 4.6; Ol has a x3.9 weaker oscillator 
strength and a xl.2 smaller wavelength. Ol is harder to 
observe due to its weaker oscillator strength, so the column 
densities considered are higher. Conversely, silicon species 
have high oscillator strengths, meaning that column den- 
sities below 10 13 cm' 2 are achievable for Sill (JV(Sin) = 

10 12.4-14.4 cm - 2) and g ilv (Ar(Sin) = 10 12.6-14.6 ^-2) 

The range for Cn is JV(Cll) = 10 132 " 15 - 2 cm -2 . Nature 
creates a convenient conspiracy whereby greater oscillator 
strengths make up for lower Type II SNe yield abundances, 
leading to similarly strong lines for our five considered ionic 
species. 

Error bars are calculated using the 

IStorrie-Lombardi et al.l (|2005l ) method (their Equation 
4), and do not include systematics such as cosmic variance. 
We also apply this error estimation method to the BSRS 
dataset, which results in much larger error bars than they 
published as a result of them including only the uncertainty 
in their summed column densities. 



current simulations owing to lack of volume. If reionization 
proceeds across such metal-enriched regions extending below 
z = 6, a dramatic decline in S7(Ol) is expected by z — 5. 

The No Field case is barely favored by the Q(Sill) data 
of BSRS, but fi(Cn) appears closer to predictions of an 
ionized IGM. The fact that observed Cn, Sill, and Ol arise 
together with similar profile shapes suggests that each of 
these ionic species may be dominant, consistent with the 
No Field assumption. If so, the lack of consistency among 
the thre e low ions may be a resu lt of metal yields differing 
those of IChieffi fc Limongil l|2004h we employ, which is one 
reason we provide alternative yields in Table [T] f2(Cll) and 
f2(Sin) show less difference between the No Field and two 
ionizing cases than f2(0 1) in Figure 1111 because the first 
two species have higher ionization potentials that are not as 
easily ionized by the addition of a field. 

Overall, the low ionization species paint a scenario in 
which there are still significant patches of neutral IGM gas 
at 2 ~ 6 along at least some LOSs. Hence in principle one 
could use such ions to trace the patchiness of reionization. 
We plan to investigate this possibility using simulations that 
more ac curately track the io nization field using radiative 
transfer jFinlator et alj|2009h . 



5.2.1 Low Ionization Species 

The resulting fi's for the five species and different fields are 
plotted in Figure 1111 We also plot the total metal density 
of the corresponding element, summed by counting the gas 
phase mass of that species in the entire simulation box at 
the given redshift (f2(species), black dash-dot line). 

With no ionization correction in the No Field case, the 
ions Cn, O I, and Sill directly trace the increase seen in 
metal density of ~ x8 from 2 = 8—^5. Note, however, that 
the ion density should be identical to the metal density if 
absorption lines of all column density are counted, yet the 
ion density is x4 — 6 less than the metal density. Where 
are the unaccounted metals? The answer is that since the 
CDDs have slope a > —2, they are in high-column sys- 
tems that are extraordinarily rare. Hence using a volume- 
weighted method as quasar LOSs to obtain a handle on a 
mass-weighted quantity such as ion density is challenging, 
more so for low ionization species, which preferentially trace 
high overdensities. 

The Q(Ol) determined by BSRS favors the No Field 
case over either ionization field, considering their sample 
from z > 6.2 quasars. This suggests that at least some of the 
z ~ 6 metal-enriched IGM is neutral. That BSRS found four 
Ol systems (all with JV(Ol) > 10 137 cm -2 ) toward SDSS 
J1148+5251 and not more than one in any of the other 8 
high-z LOSs; this is unlikely to be a statistical fluctuation 
(0.18% chance according to BSRS). The No Field case pro- 
duces a frequency of 9.2 Az _1 of such Ol lines, which is 
in agreement with the J1148 LOS and at least x20 higher 
than either of our simulated ionization fields. Hence our sim- 
ulations agree with BSRS's interpretation that the J 1148 
LOS possibly probes a long stretch of IGM with incomplete 
reionization, while the other LOSs of BSRS are consistent 
with being completely reionized. How one obtains a rela- 
tively neutral LOS along the 200 comoving Mpc spanned by 
the J1148 LOS while most LOSs are reionized is not easy 
to understand, but we cannot address this directly with our 



5.2.2 High Ionization Species 

fi(Crv) has been me asured in severa l studies at z > 4.5, 
with differing results. ISongailal <|200lh finds fi(ClV) is con- 
sistent with remaining constant from z = 2 — > 5.5 , which 
was extended by the measurements of ISimco el (120061) out to 
z = 5. 4—6.2. Conversely, the null detections bv lBecker et al.l 
(|2009h of Civ in the highest S/N data at z = 5.3 - 6.0 
supports a signi ficantly declining Q ( C iv) at z > 4.5. The 
measurements of iRvan- Weber et all (|2009l ) over the largest 
pathlength (AX = 25.1) centered at (z) = 5.76 also finds a 
similar S7(C iv) value. 

Exploring a range of wind models, OD06 explained the 
near constancy in fi(Civ) at z = 2 — 4. 5 by a decreasing 
global ionization fraction of C iv counterbalancing the in- 
creasing enrichment of the IGM. Therefore, the stability of 
fi(C iv) does not indicate that the IGM must be enriched at 
z > 6. Instead, the ionization fraction decreases with time 
for several reasons: (i) Decreasing physical densities due to 
Hubble expansion, and (ii) the increasing HM2001 ioniza- 
tion field strength from z — 5 — > 2 result in C iv tracing 
higher overdensities and less of the carbon in the diffuse 
IGM; while (iii) increasing energy input into the IGM pri- 
marily from galactic outflows at later times push metals to 
hotter temp eratures where Civ does a poor job of tracing 
carbon (see lOppenheimer et all l|2007h for a simple break- 
down of this evolution). 

OD06 predicted f2(Crv) rises by a factor of two from 
z = 6 — > 4.5 using the HM2001 field, although the new cos- 
mology we now use with a reduced as significantly curtails 
high-z SF resulting in a greater rise that agr e es be tter with 
iBecker et af] (|2009h and lRvan- Weber et all (|2009h . Figure 
1111 shows that 57(Civ) continues this growth trend between 
z — 8 — > 5, for both ionization fields we consider. For the 
HM2001 field, the C iv ionization correction falls by a factor 
of 10 (from 220 to 21) between 2 = 8^5. This is a stunning 
drop and occurs because of the ionization behavior explored 
in £14.21 the uniform background grows stronger while phys- 
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Figure 11. The evolution of the global ion mass densities, determined from simulated quasar LOSs for the three ionization cases, 
compared to the metal density summed for the corresponding atomic species calculated by summing over all the SPH particles in the 
simulation box; the difference between the two represents a global ionization correction for that ion. We sum absorbers over 2 decades of 
column density corresponding to the range the given transition is observable and unsaturated. The No Field case assumes no ionization 
correction, but falls below the total metal density because most of the absorption lies in extremely rare saturated lines, which are not 
in the included column density range when summing Q. The observations of BSRS in quasars at z > 6.2 appear to favor a partially 
ne utral IGM. Most C I V obs ervations cannot di s tingu ish between the HM2001 and Bubble models, although the rapid evolution observed 
bv lBecker et all | <2009l) and iRvan- Weber et al.l <200Sh at z > 5.3 relative to z < 4.5 supports the former. 



ical densities decline resulting in a greater ionization pa- 
rameter allowing metals to be ionized to Civ by z — 5. In 
essence, some of the same trends that cause the C IV ioniza- 
tion correction to increase from z = 5 — » 2 (OD06), cause it 
to decrease from z — 8 — > 5. 

The Bubble model predicts a slower evolution in 0(C iv) 
at z > 4.5, because this field does not evolve as much as the 
HM2001 uniform field. Metals around galaxies are ionized by 
a nearly unevolving field between z = 8 — > 4, so an increase 
of fi(C iv) is more likely to reflect an actual increase in IGM 
metallicity near galaxies. From z = 8 — > 5, f2(Civ) increases 
by xl4 while 0(C) increases by x9. As future O(Crv) ob- 
servations achieve higher redshift above 6, we predict to see 
a decline corresponding to the evolution in the global metal 
density if the Bubble field is correct. Meanwhile, the physi- 
cal and ionization conditions make z ~ 4 — 6 the ideal time 
for Civ to trace metals in the IGM for either field, as the 
ionization correction reaches a minimum at these epochs. 

O(Sirv) mirrors the behavior for O(Clv) for the two 
fields we consider, with smaller differences. The lower ion- 
ization potential of Si ill causes Si iv to trace higher overden- 
sities, which can track some of the metals already injected 
into the IGM at z = 8 with the HM2001 field. The evolu- 
tion of fi(Siiv) is nearly identical to that of 0(C iv) for the 
Bubble field for the reasons mentioned above. 



The trends in our simulated high-ionization O's ap- 
pear consistent with those that are observed. Our predicted 
fi(C iv) values with either ionization field over-estimates the 
most recent observations at z > 5.3; however we argue that 
not all metals are ionized by such a field, and that the ob- 
served value falling between the extremes of the Bubble field 
and the No Field case (i.e. no O(Crv)) is encouraging. 



5.2.3 Systematic Uncertainties in O(Crv) 

While the above agreement with 0(C iv) data for our ioniz- 
ing field models is nice, it is worth pointing out that there 
are many possible systematics in the modeling that can im- 
pact predictions. Here we point out some of the ways our 
predictions of fi(C iv) may be uncertain. 

Ionization background: Besides the possibility that 
not all metals are ionized, there is the uncertainty about the 
nature of the ionization field where it exists. The intensity 
of the ionization background at the Cm ionization poten- 
tial is unobservable and uncertain, yet is vital for the correct 
determination of the C iv ionization correction, as we stress 
in §4.2\ The carbon spectral ratio of the Bubble field makes 
O(Crv) higher than the HM2001 field; in fact, the Bubble 
field produces a nearly opti mally high value o f Q(C iv) at 
z = 6. On the other hand, iMadau fc HaardU (|2009l ) notes 



18 B. D. Oppenheimer, R. Dave, & K. Finlator 



that reprocessing of radiation between 3 and 4 Ryd by res- 
onant line absorption of the Hell Lyman series could re- 
duce by >x2 the p hotons capable of ionizing C ill to C IV. 
iBecker et alj ([2009] ) suggests this could be an explanation 
for the strong evolution they find in fi(Crv) (x4.4 increase 
between z = 5.3 — 6.0 and z = 2.0 — 4.5) if He II reionization 
proceeds over this redshift range, increasing Cm ionizing 
radiation, and resulting in a lower ionization correction for 
Civ at z < 4.5. 

Wind model: While we have investigated a plausible 
wind model for enriching the IGM, based on momentum- 
driven wind scalings, this may not be the only possible out- 
flow scenario that can reproduce observations. Varying wind 
speeds and mass loading factors in OD06 led us to some use- 
ful intuition when compared to C IV absorption data: too 
high wind speeds overheat the IGM while too low wind 
speeds do not sufficiently enrich filaments; also, too high 
mass loading suppresses SF so much that too few metals 
are produced while too low mass loading results in insuffi- 
cient metals in the IGM. But while constraints at z < 5 are 
fairly tight, at higher z the possibility of more exotic pro- 
cesses such as Pop III stars or the lack of dust may impart 
fundamental changes to outflows. Without a more detailed 
physical understanding of the physics of outflow driving, it 
is difficult to ascertain the systematic uncertainties in our 
wind model. 

Mini-halos and Population III stars: The very first 
stars and galaxies are not resolved in our simulation, as 
they likely form within mini-halos havin g virial tempera- 
tures below 10 4 K (e.g. I Abel et al.ll2002T ). If such systems 
dominate the early enrichment or io nization budg e t, our 
modeling will be incomplete. However. IWise fe Abel (|2008h 
sugge sts that they do no t dominate the ionization budget, 
and iBromm et alj ([20041 ) estimated that they increase the 
global gas metallicity by only ~ 10 -4 Zq, which is x30 less 
than the global gas metallicity at z — 6 in our simulations. 
Hence unless they have much higher yields than expected, it 
seems unlikely that metal production in mini-halos should 
be important at z — 5 — 8. 

Unresolved galaxies: For galaxies in halos with virial 
temperatures above 10 4 K, our simulation nearly resolves the 
Jeans mass at z < 10. Using the formula for Jeans length 
from §6 of OD06, we find a Jeans mass of 9 x 10 7 M at 
z — 10 assuming S = 50 and T — 10 4 K, which is nearly the 
dynamical mass limit of a resolved galaxy in our simulation. 
The flattening of the stellar mass function at M* < 10 7 Mq 
reflects the effect of the photoionization of regions containing 
forming galaxies (i.e. the filtering mass) . For this reason, our 
resolution-converged fit to the stellar mass function (Equa- 
tion [3]) suggests that we are not missing a large population 
of lower mass galaxies in the dl6n512vzw box. 

Yields: Another uncertain factor that could al- 
ter fl(Civ) is the SN e yields. We use the yields of 
IChieffi fc LimongH ([20041) , which do not change much at low- 
Z for carbon. In our simulations the SFR efficiency is so 
high in early galaxies that the average stellar metallicity in 
star particles reaches 0.1 Zq by z — 7, which minimizes the 
importance of low-Z yields. It has been suggested, however, 
that patches of low- metallicity SF persist to la te epochs (e.g. 
iTornatore et~alll2007l ; ICen fc Rique lmc 2008), which if true 
means our simulations are over-mixing metals. Integrating 
the carbon yields over the various stellar IMFs listed in Ta- 



ble [T] does indicate as much as a factor of two difference 
with the yields we u se at the highe r end. Note also that the 
yields calculated by iHirschll ([2007) for rotating stars at so- 
lar metallicity can more than double the total carbon yield 
versus non-rotating models. 

Initial mass function: Even if metal-free stars are 
unimportant, it is possible that there is an ov erall shift to - 
wards a more top-heavy IMF at high redshifts. iDavel (|2008l ) 
finds the evolution in the M*-SFR calibration could be ex- 
plained by an increasingly top-heavy IMF from z — — > 2, 
which suggests very few low-mass stars if extrapolated to 
z ~ 6. iTumlinsonl ([20061 ) presents constraints indicating a 
top-heavy IMF at early epochs that explain the frequency 
of Galactic carbon-enhanced metal poor starfl IMF varia- 
tion concerns are mitigated in large part by noting that our 
simulations match the observed UV LF, which indicate that 
they are producing roughly the correct number of massive 
stars. Hence to first order the amount of carbon produced 
will be correct. But depending on the nature and form of an 
IMF variation, it remains possible that the overall carbon 
yield per unit hi gh-mass SF cou l d stil l vary somewhat. As 
an example, the iPortinari et al. I l| 19981 ) yields using a con- 
tribution only from 30-120 Ma s t ars p roduces half as much 
carbon as the lChieffi fc Limongil ([20041 ) yields (Table Q). 

Cosmological parameters: Given the uncertainties 
on cosmological par ameters in the 5-year WMAP results 
iHinshaw etail (|2008l ). the measurement of a 8 (0.820±0.028) 
imparts the greatest uncertainty in regards to high-z SF. As 
mentioned in i]3.1l the primary difference causing the factor 
of two increase in stellar mass density at z = 6 in the w- 
series used in DFO06 over the d-series used here is as (0.90 
compared to 0.83). Similarly, we find xl.8 fewer z — 6 stars 
in a 16/i _1 Mpc 2 x 256 3 test run where we reduce as from 
0.83 to 0.75 holding all other parameters constant. Given 
the error bars on as, a factor of two difference in SF and 
chemical nucleosynthesis at z ^ 6 is possible, which will lead 
to similar uncertainty in f2(C iv) if the metals are injected 
into the IGM the same way. This is the main reason why 
the d-series results in a larger increase in fi(C iv) from z — 
6.0 — > 4.5 compared to the published w-series predictions in 
OD06. 

To summarize, a number of possibilities exist to alter 
f2(Clv) by a factor of a couple. An important constraint is 
the number of weaker lines (N(Civ) < 10 14 cm -2 ), which 
our models may be over- predicting giv e n the comparison to 
the largest sample from [Becker et all (|2009l ) that can de- 
tect such absorbers at z > 5.3; this could be a sign that 
superwinds have a smaller filling factor in the IGM than 
our simulations predict, or that more diffuse metals are less 
ionized (see H6.3.2p . Lower carbon yields, possibly from a 
top-heavy IMF, could potentially lower fi(C iv), while a dif- 
ferent ionization background is more likely also to decrease 
f2(Clv) compared to the Bubble field. The determination of 
as also can potentially alter Q(C iv) by a factor of two. The 
good agreement seen in Q(C iv) evolution should be taken 
as an encouraging but not definitive demonstration of the 
validity of our models. 



4 This itself is not a carbon yield constraint, and only an IMF 
constraint. 
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Figure 12. Predicted cumulative equivalent width distributions for absorber systems of five species at z = 8 (magenta) and z = 6 
(orange). The three ionization cases are shown, along with long dashed cyan lines for Civ and Si IV corresponding to the case where a 
uniform metallicity floor of 10 — 3 Zq is applied to the entire volume at z = 6 using the HM2001 field. The HM2001 field makes more 
strong Cn and more weak Civ systems compared to the Bubble field at z = 6. O I is the best tracer of the No Field case. Civ is ideal 
to trace a larger volume filling factor of metals, while this is not true for its sister ion Si IV, which is mostly ionized to Si V in the same 
volume. 
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5.3 Equivalent Width Distributions 

We present cumulative rest equivalent width (EW) distribu- 
tions at z = 8 and 6 for all five species in Figure fOffl This 
is intended to facilitate future observational comparisons, as 
well as present a more realistic assessment of the discrimi- 
natory power of high-z metal absorber statistics from lower 
resolution near-IR spectrographs. 

Many of the trends already discussed earlier in this sec- 
tion are reproduced here. The average weaker ionizing flux 
of the HM2001 field at z = 6 results in more strong Cll 
absorbers and fewer strong C IV absorbers compared to the 
Bubble field; the opposite is true for weak absorbers. The 
same trends are apparent in Si II and Si IV, and also at z = 8 
for Civ, although the Bubble field produces more high ion- 
ization absorbers at all equivalent widths here. Ol best dis- 
tinguishes neutral and ionized gas. 

We note that EWs, particularly at the high end, 
may be in fluenced by the phenomenon of turbulent broad- 
ening. In lOppenheimer fc Pavel | 2009) we found turbu- 
lent broadening necessary to reproduce the observed low- 
z O VI 6-parameter distribution. The purported turbulence 
occurs far below the mass scale resolved by our simu- 
lation, an d could be associated with energy injection by 
outflows l)Qppenheimer fc Pavel l2009h . Since high-z ab- 
sorbers are universally young and associated with recent out- 
flows, turbulence may result in a higher EW (for saturated 
lines), which woul d affect the high end EW distribution. 
iRauch et al.l (|200ll ) observed velocity differences in (z) ~ 2.7 
C IV profiles from lensed quasar LOS pairs at sub-kpc scales, 
which indicate more turbulence than we expect in most low- 
z Ovi absorbers. Fortunately, turbulence should not much 
affect the weak end, since turbulent broadening will leave 
EWs mostly unaltered. The other observables we have con- 
sidered such as Q(C iv) all deal with column densities, which 
should remain unaffected by turbulent broadening as long as 
lines are not saturated, fe-parameters are the observable most 
likely affected by turbulent broadening, and may explain 
some of the broad C iv pr ofiles seen by Ryan- Weber et al.l 
l|2006l ) and ISimcod l|2006l ), although iBecker et all (|2009l ) 
suggests that these may be multi-component systems with 
smaller line widths if observed at higher resolution. 

We also consider the case of a metallicity floor of 
10~ 3 Zq added everywhere in the simulation box with the 
uniform HM2001 ionization field at z — 6 for C iv (long 
dashed cyan line in the up per panel of Figure I12p. N ote that 
for a Chabrier IMF and IChieffi fc Limongil (j2004l ) yields, 
obtaining this metal floor would involve turning 0.05% 
of baryons into stars, equivalent to all the stars formed 
by z — 7.5, and distributing this uniformly among all 
baryons, which seems implausible. Hence this scenario re- 
quires a more exotic origin, possibly resulting from very 
early fast winds powered by high-yield VMSs. Such a sce- 
nario is similar to that assum e d in sem i-analytical enrich- 
ment models of IMadau et alj l|200ll) & IScannapieco et ail 
(|2002h . IMadau et all (|200ll ) calculated a complete filling fac- 
tor of metals fro m pre-galactic outflows in their most ex- 
treme cases, while IScannapieco et all l|2002l ) found as much 



5 Please contact the author for the tabular form of the predictions 
in these or any other figures. We have EW and CDD distributions 
for any other desired redshift as well. 



as 40% of the z = 6 volume is enriched to ^ 10 -3 Zq pri- 
marily by SN-driven outflows from early galaxies. 

The observational signature of such a metal floor is dra- 
matic at the weak end of the C iv equivalent width distribu- 
tion. Today's instruments should be able to integrate down 
to EW=30 mA to statistically distinguish between our stan- 
dard simulation and a metal floor: 16.9 Az _1 versus 5.7 
Az _1 . We also show this scenario in the CPP plot in the 
bottom left of Figure \TU\ (long dash ed cyan line). The high- 
est S/N observations yet at z > 5.5 (I Becker et al.|[2009l ) find 
no weak components, indicating that a significant fraction 
of the IGM is not enriched far above 10 -3 Zq. Pata is cur- 
rently inconclusive as to the extent of volume enrichment 
at these low levels, but sitting on the brightest quasar at 
2 > 6 for several nights on an 8-10 meter telescope should 
produce the necessary S/N of ~ 40 — 50 to adequately test 
this scenario. C iv is by far the best species for testing such 
a scenario, because it is the most effective tracer of metals 
around mean overdensity at z — 6; Siiv is ionized to Siv at 
the same overdensities and the metallicity floor does affect 
the weak absorber frequency. 

Finally, we note a contrast in the C iv distributions at 
high-z to the z = — 0.5 EW distribution for O vi. The faint 
end turnover of the latter sho ws very little sensitivity to th e 
volume filling factor of metals (jOppenheimer fc Pave|[2009l ). 
Part of the difference is that Ovi is not tracing the lowest 
overdensities where metals reside at low-z; the C iv contours 
overlap the metals at the lowest overdensities as shown in 
Figure [7] This is precisely the reason why C iv holds so much 
potential for mapping the true extent of high-z metals. 



5.4 Aligned Absorber Ratios 

Another observationally-accessible statistic is aligned ab- 
sorber fractions. We consider absorbers of different ions to be 
aligned if both components are stronger than 50 mA and are 
within 15 kms -1 of each other (5v ^ 15 kms -1 ). We exam- 
ine the No Field first in Figure 1131 The alignment fractions 
are listed in the corners, where the percentage corresponds 
to the likelihood of finding an absorber of the species on 
the opposite axis aligned with the absorber species of the 
near axis. For the neutral case, the alignment percentages 
are above 95% in every case. This lends support to the idea 
that the low ionization species observed by BSRS, which ap- 
pear aligned in every case for O I observed above z = 5.8, 
are tracing metals where there is no ionizing flux above the 
Lyman limit. These systems occur far too frequently to be 
arising within the ISM of intervening galaxies according to 
BSRS. This could be evidence that patchy reionization is 
still occurring at these redshifts, or alternatively that dense 
cloudlets entrained in winds are self-shielded from the ioniz- 
ing background. The latter scenario would predict that there 
be high ionization species in the same systems, possibly off- 
set, indicating multi-phase outflows a s observed locally (e.g. 
IStrickland et alj|2002t iMartinl l2005bl) a nd seen in galactic- 
scale simulations (e.g. lFuiita et al.l l2008). Therefore C iv and 
Si IV could be very useful additional constraints. 

The high alignment fractions are helped by the fact that 
these three species have similar EWs for typical Type II 
SNe yields. Aligned absorbers fall along the linear relations 
in each subplot of Figure 1131 where th e ratios in column 
densities equal the ratios in Type II SNe lChieffi fc Limongil 
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Figure 13. Aligned absorbers of different ions are denned as having both components above 50 mA and a velocity difference less than 
15 kms -1 . Here we present the No Field case, which at z = 6 shows ab sorbers fall along a linear relation with some scatter; the lines 
show the expected relation if Type II SNc yields (Chicffi & Limongi 2004, Z = 10~ 3 ) with no ionization corrections assumed. Alignment 
fractions are listed in the corners, where the percentage corresponds to the likelihood of finding an absorber of the species on the opposite 
axis aligned with the absorber species of the near axis. The percentages exceed 95% in every case, and remain similarly high when the 
velocity tolerance is reduced to 8 kms -1 , because they are tracing metals in the same neutral gas. 



(|2004h yields. Although we use a component separation 
criterion of Sv 15 kms" 1 , the alignment statistics are 
nearly identical when Sv = 8 kms" 1 and Sv = 32 kms" 1 . 
This is generally the case in all alignment statistics consid- 
ered here, which indicates that the species are tracing the 
same underlying gas. This contrasts with, for instance, low- 
z O VI — Hi alignment fraction s which are markedly l arger 
at higher Sv (|Tripp et al.ll2008l ; iThom fc Chenll2008al ). be- 
cause O VI traces physically d istinct gas from the Lya-forest 
|Oppenheimer fc Davell2009] ). 

The alignment between low and high ionization species 
in Figure [T4l indicate many fewer aligned absorbers. Cn is 
aligned with C IV 45-62% of the time, while C IV is found 
aligned 11-23% of the time with C II. Finding C II where C IV 
exists should occur much less frequently than the other way 
around according to the visualizations in Figures [S] and [5] 
the small filling factor of C n is more often coincident with 
C iv, while C iv can also trace the more extended metals. 
This is an ideal ratio to measure the shape of the ionization 
background, since abundance ratios are irrelevant and the 
carbon spectral ratio that is measured covers a large ioniza- 
tion range from 11.3 to 47.9 eV. There exist stronger Civ 
components aligned with C n for the Bubble field relative to 
the HM2001, due to t he greater hardness of the former field. 
In the one component iRvan- Weber et all (|2006l ) could mea- 
sure both at z = 5.7239 in SDSS J1030+0524, they found 
no aligned C II. A sample size of at least a dozen strong 
C IV components exploring aligned C II should begin to build 
enough statistics to constrain the shape and intensity of the 
ionizing background. 

Si II and Si iv alignment provides a similar opportunity 
as the carbon species to measure the shape of the back- 
ground, and the alignment fraction is even higher than for 
carbon due to the lower ionization potential of Si ill (33.5 
eV) compared to C ill. BSRS actually does find aligned Si iv 
in three of their components lying close together in the 
2 = 5.3364 system of SDSS J0231-0728. This system in- 
dicates possible non-uniform ionization, because the very 
aligned profiles of Cn, Ol, and Sill in components about 
100 kms" 1 away from the aligned Sill — Siiv is a strong in- 
dicator of the No Field case, while Si iv requires an ionizing 
field. BSRS mentions that aligned C iv is also present ac- 



cording to lBarth et al.l |2003t ). but the low resolution of that 
observation makes it difficult to determine the component 
alignment; we expect it to be aligned to the components ex- 
hibiting the Si iv. The presence of low and high ions nearby 
suggests a galactic ISM as might be seen in a damped Lya 
system, as forwarded by BSRS. 

C n is usually found with O I, but the reverse is not 
true if there is an ionization field above the Lyman limit. 
The bluest component of the z — 5.3364 BSRS system men- 
tioned above is consistent with a non-detection of O I with 
the observed C n and Si II, supporting some sort of ionization 
field. Finding the fraction of C n absorbers with aligned O I 
is one of the most sensitive indicators of the three different 
ionizing cases we explore. 

Finally, we do not show the incidence of O I aligned with 
C iv, because this case rarely occurs under an ionizing field 
due to the large difference in ion ization potentials. T his is 
apparent in SDSS J1148+5251 as lBecker et~aTI (|2009T i finds 
no C iv despite the 4 O I absorbers (BSRS) . 

Overall, the alignment fractions offer another way to 
constrain the nature of the ionization field. Current data 
are in line with expectations from our models, though this 
is not highly constraining. Larger samples will elucidate the 
nature of individual absorbers, and possibly allow one to 
measure the carbon spectral ratio which could give clues as 
to the nature of ionizing sources. 



5.5 System Profiles 

Finally in this section, we consider the appearance of sys- 
tems in absorption line spectra. As an illustration, we choose 
a series of LOSs through the region indicated in Figure 1151 
which is the same region shown evolving in Figure|5]centered 
on the second most massive galaxy in the simulation box. We 
choose two LOSs that go exactly through the center of this 
M* — 10 9 ' 1 M© galaxy and show the system profiles for our 
five ions and our three fields in Figure 1161 The a;-direction 
spectrum is more along the direction of the local filament, 
while the y-direction is nearly along the axis of the bipo- 
lar outflow. Both show extremely strong systems, but the y 
spectrum is remarkably wide (800 kms" 1 ) owing to pecu- 
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Figure 14. Aligned absorber ratios and percentages displayed analogously to Figure IT3l for the HM2001 and Bubble ionization fields. 
Symbols indicate components where both absorber species are > 50 mA and Sv < 15 kms -1 , while percentages correspond to the 
likelihood of finding an absorber over 50 mA of the species on the opposite axis aligned with the > 50 mA! absorber species of the 
near axis. Aligned ratios can help distinguish the shape and magnitude of the ionization field, especially in the case of C II — Civ. The 
alignment fraction of C II with O I is one of the best ways to distinguish the three ionization cases, although these are rare if the IGM is 
ionized. 
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liar velocities imprinted by the enriching outflows. Compar- 
ing the physical parameters of the gas in the bottom three 
panels without peculiar velocities (cyan) and with peculiar 
velocities (red) shows that the velocity structure dominates 
the system profiles of both LOSs. This galaxy has a maxi- 
mum v w i n d just over 500 kms -1 , therefore Sv ~ 400 kms -1 
on each side of the y spectrum is sensible as winds are 
launched with an intentional spread of velocities, plus the 
LOS is slightly offset from the outflowing axis. 

None of the observed systems in the data (or in our 70 
LOSs) have profiles as strong and as wide as the systems 
in Figure 1161 because not surprisingly lines of sight directly 
though the centers of galaxies are rare. If we assume similar 
profiles arise within 3 kpc of the 10 galaxies in our sim- 
ulation box over M„ = 5 x 10 8 Mq, then the frequency of 
intersecting such a system occurs once in Az = 2000! A more 
relevant consideration for the spectra in Figure [16] is to take 
the blue-shifted half and model the outflows in the spectra 
of LBGs at lower redshift (z ~ 2 — 3), where such galactic 
superwinds have been obs erved to be common and observa- 
tions are more obtainable dPettini et al.ll200ll ; IShaplev et al.l 
2003); this is beyond the scope of this paper. We note that 
the faster moving materials absorb more in the lower ion- 
ization species with the weaker HM2001 background. 

More typical examples of quasar absorption line systems 
are shown in Figure [17] where two LOSs have impact param- 
eters, 6, of 30 kpc to the same galaxy. Velocity profiles of ab- 
sorbing ions are still dominated by peculiar velocities, and 
the LOS closest to the outflow axis (y direction) produces 



Figure 15. The 4 X 4/i _1 Mpc X 25 km s~ 1 region shown in Figure 
[9] at 2 = 6 with intersecting LOSs, of which we show the absorp- 
tion line profiles in Figures 1161 and 1171 This snapshot spans 950 
physical kpc, and the Hubble parameter is 651 kms -1 Mpc — 1 . 
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Figure 16. The LOSs in Figure Il5l intersecting the center of the Mr = 10 9 ' 1 Mq galaxy. The y LOS (right) shows more extended 
absorption line profiles (Sv ~ 800 kms^ 1 ), because it is along the bipolar outflow axis, while the x LOS (left) follows more closely the 
local filament. No systems in the data or in our 70 simulated LOSs match such profiles, because intersecting such regions is extremely rare 
when considering quasar LOSs, which probe a volume-weighted sample. The blue-shifted side of these profiles may be more appropriate 
when considering the absorption of the outflows and galactic ISM in the spectra of LBGs and high-z galaxies. Lower ionization species 
show broader profiles in the HM2001 case versus the Bubble case. 



a wider system. This system is perhaps most qualitatively 
similar to the two wide systems observed in low ionization 
species by BSRS where Cn, Sin, and in one case Ol show 
a total velocity spread of Sv ~ 250 kms - . This supports 
the idea of neutral/low-ionization gas entrained in outflows, 
self-shielded from the ionization background; this probably 
makes more sense for the z = 5.3364 system than patchy 
reionization extendin g to this redsh ift. A similar velocity 
spread is observed by ISimcoel (|2006l ) in his Civ z = 5.829 
system. 

Although an alignment within 30 kpes occurs once in 
Az = 20 (if we count the ten most massive galaxies in the 
dl6n256vzw box), we will argue in i]6.3l that such systems are 
qualitatively similar to systems observed around lower mass 
galaxies at lower impact parameters. The frequency of in- 
tersecting such qualitatively similar systems as those shown 
in Figure [T7J is Az ~ 2.8 if we consider all galaxies M* > 

10 7 Mq, and d ga i < 30y ygfr. Two such systems should 

ha ve been inters e cted b y n ow given the combined pat hlength 
of iBecker etafl (|2009D fc iRvan- Weber et all (|2009D . which 
compares well to three strong C IV systems being observed. 



Three differences when scaling to lower stellar mass are (i) Sv 
of systems are smaller as the outflow velocities are smaller, 
(ii) metallicities will be lower due to the mass-metallicity 
relation, and (iii) the mass loading factors should be higher 
at lower mass, which counterbalances (ii). 

In Figure[l8]we show how median Sv varies as a function 
of column density in our 70 LOS z — 6 sample. Here Sv is 
defined as the difference between the central velocities of the 
two furthest components in a system (Sv = when there is 
a one component system ) . One a boundaries (dashed lines) 
show the range of velocity spreads in each case. This can be 
compared to a sample of the same spectra without peculiar 
velocities, shown as long dashed orange lines. While this is 
not strictly a fair comparison as systems without peculiar 
velocities will not necessarily fall into the same column den- 
sity bin, this nevertheless clearly illustrates that the stronger 
observed systems are dominated by peculiar velocities. At 
high column densities, greater than 50% of a system's width 
owes to the peculiar velocities from outflows. With future 
samples, such a plot can be used to infer conclusions about 
the imprints of outflows on high-z metal-line profiles. We do 
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Figure 17. The two LOSs illustrated in Figure IT51 with impact parameters of 30 kpc. As in Figure [T6l peculiar velocities dominate the 
system width, and the y LOS (right) exhibits a more extended profile as this is close to the axis of the bipolar outflow. These syst ems are 
more similar to those observed in quasar LOSs, and could be analogous to some of the broader systems observed by BSRS and ISimcoei 
j2006h . 



not make a direct comparison now, because of the paucity 
of observed systems and differences in both S/N and reso- 
lution, which may affect such a comparison. 

The main points here are that (i) we produce systems 
qualitatively similar to those observed in the data at a fre- 
quency that seems reasonable; and (ii) the velocity widths of 
the strongest observed systems are likely to be dominated by 
peculiar velocities associated with outflowing gas. The lat- 
ter suggests that plots of Sv vs. column density may provide 
constraints on outflow kinematics from high-z galaxies. 



6 THE PHYSICAL ENVIRONMENT OF 
ABSORBERS 

As our simulations form galaxies and enrich the IGM to- 
gether, we are able explore the relationship between galaxies 
and metal absorbers, along with their physical and environ- 
mental properties, specexbin tracks for each absorber the 
physical properties of density, temperature, and metallicity; 
the environmental properties of the distance and stellar mass 
of the nearest galaxy; plus the last time an SPH particle was 
launched in a wind (i.e. its "age"). We also track an SPH 
particle's SFR if it is in a galaxy, or « w ind at which it was 



launched if in an outflow. Unless otherwise mentioned, we 
use the Bubble model, because (i) the ionizing flux from the 
local galaxy should dominate at the distance of the metals 
given the reasonable assumptions for the transmission above 
the Lyman limit, and (ii) this field produces a shallower col- 
umn density distribution, which seems to be supported by 
the data. Note, however, that the No Field approximation 
may best explain some absorbers, most notably those in the 
J1148+5251 LOS. This may help remedy the overestimated 
Bubble Q(C iv) if more metals are at lower ionization states 
as we suggest in H6.3.2I We also address the physical environ- 
ment of the HM2001 field, which may be most appropriate 
at z < 6, in flolTTI 

We concentrate on C IV absorbers at z = 6, examin- 
ing first their physical state and then their evolutionary 
status, and bring them together by exploring the high-z 
galaxy- absorber connection. We consider an absorber's con- 
nection to the stellar mass of its originating galaxy, which 
in our models is a good proxy for the more observationally- 
accessible SFR. This leads us to a suggestion of possibilities 
of why we may produce too many weak absorbers at high-z. 
We then switch to lower ionization species represented by 
C n, and lastly extend our analysis back to z = 8. 
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Figure 18. System width (8v) versus column density for C IV 
HM2001 and Bubble fields (top and middle), and the O I No Field 
case (bottom). Solid lines indicate the median and dashed lines 
show 1 o~ dispersions. Sv is defined as the velocity difference be- 
tween the central velocities of the two furthest components in 
a system. The orange long-dashed lines show 5v if we make our 
LOSs without peculiar velocities (i.e. the spatial distribution con- 
tributes to the system width). The majority of the system width 
is due to peculiar velocities, where outflowing metals broaden the 
profile. 



6.1 Physical Conditions 

Figure [19] shows four absorber properties plotted against 
each other for z = 6 C IV absorbers in the Bubble field. In the 
upper left panel are the same type of logarithmic metallicity 
contours considered in Figure [7] C IV has the smallest global 
ionization correction at z ~ 5 — 6, because the ionization 
fractions of this species are ideally aligned with the range 
of densities at which metals reside at this redshift (8 ~ — 



100). Histograms in the bottom left panel show that stronger 
absorbers trace higher overdensities. 

Temperatures indicate predominantly photo-ionized 
C IV. Outflows typically heat gas to 10 5 " 6 K around galax- 
ies, but the gas quickly cools to warm IGM temperatures 
due to efficient cooling in this range. This appears to con- 
tradict our previous work in OD06, where we claimed a sig- 
nificant fraction of z ~ 6 C IV (43%) is collisionally ion- 
ized, but the difference arises because in OD06 we quoted 
a mass- weighted fraction (from particles) , while here we are 
computing a volume- weighted quantity (from spectra). For 
8 < 200, our present models agree well with OD06, com- 
paring Figure 9 in OD06 with the upper left of Figure 1191 
Regions with 8 > 200 are mostly forming stars at z = 6, and 
the frequency of metal systems where N(C iv) 10 13 cm -2 
includes SPH particles above the SF threshold is 0.4A2 - . 
Of the 54 components with AT(Civ) > 10 14 cm" 2 , only 12 
have any traces of SF. Ironically, only one of these 12 lines 
has its line center density at 8 > 200; even when LOSs probe 
galactic environments, the C IV absorbing gas is primarily in- 
tergalactuQ. Hence our models clearly predict that the vast 
majority of typical C IV absorbers are primarily intergalactic 
at z — 6, and the rest mostly trace halo outskirts. 

The average carbon metallicities at iV(Civ) = 10 13 
and 10 140 cm -2 are 0.05 and 0.10 Zq respectively. This is 
about ~ x 15 — 20 the IGM metallicities at the correspond- 
ing overdensities, using the simulation-averaged metallicity- 
density relation in Figure [3] Hence C IV absorbers are trac- 
ing a clumpy, inhomogeneous distribution of metals with a 
very low volume filling factor, < 0.1%. This is generic in 
our models , and is even seen at the p resent day for O VI 
absorbers (jOppenheimer fc~D avc 2009). The metallicity in- 
homogeneity is even more extreme at z = 8: 0.06 Zq for 
iV(Clv) = io 13 ' 5 - 14 cm -2 , while the simulation-averaged 
metallicity is 9 x 10 -4 Zq (a factor of ~ 70). This shows that 
metal-line absorbers trace an increasingly clumpy distribu- 
tion of metals at early times going back from z — 6 — > 8, a 
trend that extends to low- 2; IGM metals as traced by O VI 
jQppenheimer fc Davell2009l ). 

In summary, C IV absorbers seen at z ~ 6 are mostly 
expected to be true intergalactic absorbers, not gas within 
galaxies or galactic halos. Galaxy outflows reach a larger 
comoving volume at early times (OD08), but in actuality 
they are enriching a small fraction of the volume in a highly 
patchy fashion. 



6.2 Origin of IGM Metals 

We now consider how an absorber's age relates to its envi- 
ronment, as a way of characterizing in which regions IGM 
metals originate. An age-density anti-correlation exists in 
z — 6 C IV, with a significant spread, as shown in the bot- 
tom left panel of Figure [19] This is the same trend we see in 



6 The one absorber tracing primarily galactic gas has <5 ~ 1000 
and JV(Clv) = 10 17 ' 5 cm -2 ! This may be an unphysical amount 
of Civ, but this one absorber comprises 88% of the total f2(C iv) 
along all LOSs. This example illustrates the reason for the 
discrepancy between the volume-derived O(Civ) of jV(Crv) = 
igl3. 0—15.0 cm — 2 absorbers versus the mass-derived f2(Clv), as 



discussed in i|5.2l 
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Figure 19. Physical conditions of Cw. C IV absorbers at z = 6 with the variable Bubble field are plotted in planes of density, temperature, 
absorber age, and distance to the nearest galaxy. Logarithmic contours at 0.5 dex steps in the p — T phase space correspond to the 
metallicity- weighted density (darker contours are greater densities). In the d ga ;-age plane, an effective outflow velocity can be calculated 
by dividing the two quantities; we show a solid gray line corresponding to v e g = 100 kms -1 , with the 30 and 300 kms -1 indicated 
by the dashed gray lines to the left and right. Histograms along the side show the summed SiV(ClV), with each color corresponding 
to the column density range indicated by the key in the upper right panel; all absorbers where iV(Clv) > 10 15 cm -2 are counted as 
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arise from higher overdensities and are generally younger. Most absorbers are 30 — 300 Myr old, lying 5 — 50 kpc from their parent galaxy. 
The average v c ff ^ 100 kms ^ far from galaxies, which is half as much as the typical ejection speed ^wind- 



low- z Ovi (| Oppenheimer fc Dave|[2009l ) where stronger ab- 
sorbers trace higher densities that have been enriched more 
recently. In the case of C IV, absorbers over 10 14 cm -2 are 
more likely to trace metals ejected on average 100 Myr ago, 
while weaker absorbers trace metals 100-300 Myr ago from 
earlier galaxies when they were at lower mass and metallic- 
ity. 

The age-density relation can be understood by consid- 
ering the age-galaxy distance (d ga i, shown in physical kpc) 
relation in the bottom right panel of Figure 1191 Older ab- 
sorbers have had time to travel greater distances reaching 



lower overdensities. Taking d ga i divided by age, we can plot 
the average velocity these winds need to travel to get as far 
as they do, which we call the effective wind speed, v e s- The 
gray solid line shows v c h — 100 kms -1 and the dashed lines 
on the left and right sides are 30 and 300 kms -1 respec- 
tively. Comparing this to the average wind speed at launch 
fwind = 200-250 kms -1 (cf. Figure 0, we see that outflow- 
ing gas typically slows 20 — 50% by the time it reaches the 
IGM. As argued in OD08, this slowing is primarily hydro- 
dynamical as outflows run into a dense IGM, though gravity 
plays a non-trivial role. 
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The good age-d ga i correlation indicates that these ab- 
sorbers trace metals leaving galaxies for the first time. This 
also is expected from the evolutionary picture of feedback 
from z — 6 — > presented in OD08, where we found that 
metals more often than not recycle (i.e. return to a galaxy at 
later time and are shot out in another wind/form into stars). 
Metals recycle on a timescale of 1-3 Gyr reaching a median 
farthest distance of ~ 60 — 100 proper kpc from their parent 
galaxies. C IV is tracing metals on their first journey into the 
IGM, because the ages and distances in Figure[l9]are nearly 
all less than the recycling times and turn-around distances. 
These absorbers primarily reside i n the IGM and not in ha- 
los, b ecause the halo radius (from [Navarro . Frenk. fc White] 
119971 ) is < 5 physical kpc for all z = 6 galaxies in our simu- 
lation; nearly all absorbers are at distances above this value 
in Figure 1191 Hence C IV absorbers are tracing IGM metals 
on their first journey outward from galaxies at z = 6. 

As another way of considering absorber origin, we plot 
the ratio of the effective velocity to the launch velocity 
(ueff/vwind) versus the distance from a galaxy (left panel) 
and the launch velocity (right panel) in Figure [20] This ratio 
should not exceed one, given that wind particles are slowed 
as they travel outward. A few absorbers at greater than one 
are mostly due to the rare instances when the neighboring 
galaxy is not the originating galaxy for a wind. The clearest 
trend is that Weff/^wind shows a general decline with d ga i. For 
absorbers with N(C iv) = 10 1 3 5-14 5 cm" 2 , the s treng th of 
most compone nts observed bv lRvan- Weber et al.l (|2006f ) and 
ISimcoei l|2006h . the median ratio is nearly unity at 3-5 kpc, 
0.8 at 10 kpc, and 0. 5 by 50 kpc. Note that our ou tflow pre- 
scription (following ISpringel fc Hernquistl [2003a ) turns off 
hydrodynamic forces for a time that does not exceed that 
required to cover 4 physical kpc (at z = 6) at v w ind- This 
shows that, as soon as hydro forces are turned on, outflows 
begin slowing down as they travel further from their parent 
galaxy. 

Several other physically interesting trends are also ev- 
ident in the right panel of Figure 1201 First, stronger ab- 
sorbers arise from more massive galaxies with greater w w i n d; 
more massive galaxies enrich the IGM with higher metallic- 
ity outflows due to the mass-metallicity relationship at high 
redshift (Z oc M° a ]; DFO06). Second, winds from massive 
galaxies slow down quicker as indicated by lower v c a/v w i n d 
at higher i> w i n d; this arises because massive galaxies tend to 
be surrounded by denser gas (OD08). 

6.3 The Galaxy- Absorber Connection 

While we have explored numerous relations between ab- 
sorber physical and environmental parameters at z = 6, 
there is as of yet no over-riding correlation to explain C IV 
absorber strengths, as each parameter exhibits considerable 
scatter versus iV(Civ). Here we argue that the fundamen- 
tal relation governing metal absorbers is the stellar mass, 
M*, of the neighboring galaxy versus the distance to that 
galaxy, d ga i , and most other correlations are merely a conse- 
quence. The neighboring galaxy is defined as the galaxy with 
the greatest impinging flux on an absorber (i.e. the nearest 
galaxy weighted by the inverse square law). Figure [2T1 shows 
the M* — d ga i relation. As is evident, the relation between 
dgai and N (C iv) at a given M„ shows a strong correlation. 
We plot gray dotted lines corresponding to the rela- 



tion dgai oc V-M* representing the constant flux levels re- 
sulting from the Bubble field due to the inverse square law. 
If all physical conditions are equal, absorbers along these 
lines should be identical. Indeed, this seems to hold reason- 
ably well. Hence the key deterministic relation for absorber 
strength is distance from galaxy of a given stellar mass, more 
so than any of the physical parameters explored in the pre- 
vious sections. The paucity of absorbers below the bottom- 
most gray dotted line is only a selection effect poorly sam- 
pling these regions in a volume- averaged measurement; very 
strong absorbers exist here if LOSs are fortunate enough to 
sample them (e.g. the extreme systems displayed in Figure 

Absorbers are closely related to galaxies at high-z 
moreso than any other time because (i) metals are on their 
first journey into the IGM, (ii) they are still in close proxim- 
ity to their parent galaxies, and (iii) galaxy evolution is rela- 
tively simple to understand at high-z as galaxies are mostly 
evolving in isolation. We use these assumptions to construct 
the solid gray "evolutionary tracks" that metals take into 
the IGM, assuming constant v c g — 100 kms" 1 and an ex- 
ponentially increasing M, with time, which approximately 
falls out of the integration of Equation [4] Although both 
of these assumptions are not exactly true, this treatment 
reveals some important behaviors. 

First, "isochrones" travel parallel to the x-axis of Fig- 
ure [5T] by construction if « cn = is assumed to be independent 
of M* , meaning the absorbers of the same age are found at 
a similar distance. This of course an oversimplification, be- 
cause v a g slows down as time passes and winds with higher 
"wind slows down faster. However, this explains why there 
are very few weak absorbers (iV(Civ) < 10 13 cm -2 ) around 
M t > 10 s Mq galaxies; metals simply have not had time 
to travel far enough to the low-density regions where weak 
Civ absorbers would arise. 

Second, tracing the evolutionary tracks back to the bot- 
tom reveals the mass of the launch galaxy for each metal 
absorber. DFO06 plotted the mass-metallicity relationship 
of z = 6 and z — 8 galaxies, finding that gas metallicity 
increases by ~ 0.3 dex per decade of M*. Although the over- 
riding trend for absorber strength depends on the ionizing 
flux in the Bubble field, a more subtle trend (hardly visible) 
is that absorbers grow stronger rightward along the same 
dotted gray line primarily as a result of the mass-metallicity 
relationship. However, this effect is about half as strong as 
expected (~ 0.15 dex), and should be counter-balanced by 
increased mass loading factors of smaller galaxies. Also, con- 
sider that v e g is not really constant, which limits us from 
over-interpreting this plot. 

The only candidate observed high-z galaxy-C iv ab- 
sorber connection is along the J1030+0524 LOS, where 
iBecker et al.l (|2009T ) notes the observed systems at z = 5.7 — 
5.9 may be rel a ted to an excess of i-band dropouts observed 
by iKim et alj (|2008T ) along this quasar LOS. In our mod- 
els, most absorbers between iV(Clv) = io 13 0-15 cm -2 
are associated with galaxies where M* = io 7 0-8 - 5 Mq, 
which have corresponding 1350 A AB magnitudes rang- 
ing between —15.4 and —18.9 (see the luminosity range 
bounded by green dotted lines in Figure [2}, and SFRs from 
0.05 - 1 M yr _1 . The closest of the J1030 absorbers lies 
about 50 arcseconds away from a galaxy, i.e. an impact pa- 
rameter of 300 kpc, which is too far for our winds to reach 
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Figure 20. The ratio of a wind particle's effective velocity (d ga i divided by age) and launch velocity t) w ind is plotted against d ga i 
(left) and f w i n( j (right). This ratio generally declines the farther it travels from a galaxy as hydrodynamic and gravitational forces slow 
outflowing metals. Stronger absorbers arise from faster winds launched from more massive galaxies, but these winds are more likely to 
slow down faster as indicated by a declining f e ff/^wind at higher f w i n d- The color key for absorber strength in Figure [T9l applies here. 




log[M. (M Q )] 

Figure 21. Galaxy-absorber connection plot: d ga i vs. M* elu- 
cidates the fundamental and close connection between galaxies 
and absorbers at high-,2. It reveals a strong correlation between 
./V(Civ) and d ga i for a given mass galaxy. The ionizing conditions 
are identical in the Bubble field along the gray dotted lines, ex- 
plaining much of this correlation. Evolutionary tracks (solid gray 
lines for winds emanating from a Af„ = 10 6 , 10 7 , & 10 s Mq 
galaxy) tracing an absorber back to its parent galaxy assuming 
a constant v e ff = 100 kms - 1 and secular galaxy growth. Cur- 
rently observable z = 6 C IV absorbers are launched from galaxies 
between M, = lo 6 - 5 " 8 M Q , but are observed around galaxies 
with Af* = 10 7 0— 8 ' 5 Mq. The color key for absorber strength in 
Figure [T9l applies here. 



(e.g. Veff > 1000 kms -1 for 300 Myr); plus the magnitudes of 
these galaxies appear to be brighter than any of the galaxies 
in the dl6n256vzw simulation box. This may be an extended 
overdense region of space where smaller, currently unobserv- 
able galaxies are primarily responsible for the enrichment. 



6.3.1 HM2001 Background Case 

The uniform HM2001 background case may be more appro- 
priate below z — 6, and actually better explains the observed 
frequency of strong C IV absorbers at (z) = 5.75, although it 
is more discrepant in over-estimating weaker absorbers. We 
show the galaxy-absorber connection plot in Figure [22] for 
this field and note the following differences. First, the dis- 
tance and mass of the neighboring galaxy matters less for 
the characteristics of these absorbers; the dotted lines do 
not have any meaning, but are shown for relative contrasts 
with Figure 1211 As this field is weaker near galaxies, there 
are fewer strong absorbers because carbon is in lower ioniza- 
tion states. For the same reason, there are more weak C IV 
absorbers at greater distances from galaxies, because these 
metals are otherwise ionized to C v by the Bubble field. 

A potentially accessible distinguishing characteristic 
of the two fields involves components with iV(C iv) > 
10 140 cm -2 , which are only found around M* ~ 10 9 M Q 
galaxies in the uniform case. These absorbers arise in the 
situations like those illustrated in Figure [TBI where a bipo- 
lar outflow copiously enriches the low density IGM outside 
of filaments where the HM2001 C IV ionization correction is 
lowest. The strongest absorbers actually arise from overden- 
sities less than 10, which is the opposite of the behavior in 
the Bubble model. 
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Figure 22. The galaxy-absorber plot for the uniform HM2001 
model, which is analogous to Figure l2l1 for te Bubble model. Here 
the gray dotted lines do not have any significance, but are shown 
for reference. Fewer strong absorbers exist near galaxies and more 
weak absorbers exist far from galaxies due to the weaker ioniz- 
ing field in such environments. One important difference with 
the Bubble field is that strong absorbers are only associated with 
M* ~ 10 9 Mq galaxies, and these absorbers trace enriched low- 
density regions further from galaxies. The color key for absorber 
strength in Figure [T9l applies here. 

6.3.2 Where Are The Weak z — 6 Civ Absorbers? 

Both our backgrounds produce a significant amount of weak 
absorbers (JV(Clv) = 10 13 ~ 14 cm" 2 ), contributing non- 
trivially to fi(Clv) (30% in the Bubble case and 70% in 
the HM2001 case). The two observational z = 6Civ ODD 
constraints shown in Figure[l0] (bottom right) suggest a pos- 
sibly shallower distribution than either of our models pro- 
duce, with the C IV mass density coming almo st exclusively 
from strong systems. The iBecker et al.l (|2009l ) data point is 
a conservative upper limit, and already provides a powerful 
constraint that we interpret as indicating a significant vol- 
ume of the IGM is not enriched far above Z = 10 -3 Z@. 
However, if we consider that they observe no C IV absorbers 
between iV(Clv) = io 13 - 5-140 cm~ 2 and assume their com- 
pleteness is 75% over this column density range (this seems 
reasonably conservative using their Figure 3), the difference 
with our models is even more stark. We predict 3.3 and 3.9 
such absorbers should be detected at z — 5.65 for the Bubbl e 
and HM2001 fields respectively over the lBecker et al.l (|2009l ) 
adjusted pathlength of Az — 1.65. 

A possible explanation is that the small galaxies in our 
simulation do not produce C IV absorbers for some reason. 
The evolutionary tracks in Figure[2I]can be used as a stellar 
mass cut for which galaxies produce C IV systems. For ex- 
ample, if the originating M* must be ^ 10 7 Mq, then only 
absorbers right of the middle evolutionary track contribute. 
Systems below N(C iv) = 10 14 cm -2 are more than halved, 
while leaving the stronger end unaffected if applying this 
cut. The power law slope now approaches a = — 1 (i.e. equal 
numbers of absorbers in low and high density CDD bins). 

Why would low-mass haloes not produce C IV ab- 



sorbers? A possible explanation is that the highly mass- 
loaded outflows from small galaxies self-shield metals lead- 
ing to them being observed in lower ionization states. Thus 
the over-excess of C IV absorbers and fi(C iv) in either ion- 
ization background may be counter-balanced by the BSRS 
low-ionization systems, which are underestimated when ap- 
plying these backgrounds to all metals. 

A second possibility for us overestimating weak C IV 
absorbers is that our outflows are traveling too far into the 
IGM. For example, there exist few very weak, N(C iv) < 
10 13 cm -2 , absorbers around galaxies M„ > 10 s Mq in Fig- 
ure [2H because winds have not had time to enrich regions 
at dgai > 100 kpc. If wind speeds are even lower, metals 
will populate less of the M* — d ga i phase space at large d ga i, 
where they otherwise would produce weaker C iv absorbers. 

Our discussion in this section on the galaxy-absorber 
connection may be currently inaccessible to observations 
and somewhat speculative, as only a handful of absorbers at 
z ~ 6 are known thus far and identifying the nearest galaxy 
may have to await a next-generation facility. Our purpose is 
to emphasize the M* — d ga i relation as the key one for un- 
derstanding the origin and subsequent evolution of high- z 
metal-line absorbers. As these early galaxies push a major- 
ity of their metals into the IGM, the metals are still close 
enough that the high-z galaxy-absorber connection tightly 
establishes the absorber properties. This is a fundamental 
prediction of our model, and if proved correct, offers a novel 
way to probe feedback processes in the reionization epoch. 

6.4 Low Ionization Species 

We now examine the physical conditions of C II absorbers, 
as a representative of lower ionization species. In H4.2I we 
showed C II traces higher overdensities closer to galaxies; 
therefore we might expect this species to trace metal- 
enriched outflows at an earlier stage. This turns out to be 
only partially true. In Figure [23] we plot C II absorbers in 
p — T phase space demonstrating that C II absorbers do not 
trace the same extent of this phase space as C iv. Our vi- 
sualization in Figure shows C II tracing less volume than 
C iv. Another difference is there is a stronger correlation 
with density for C II, much more so than for C iv. The rea- 
son for this is that the C II ionization fraction falls off rapidly 
toward lower overdensities as shown in Figure [7] while the 
C iv ionization fraction peaks at overdensities correspond- 
ing to intergalactic metals. C II absorption strength depends 
much more on the physical state (density) than C iv, which 
depends more on its environmental state (proximity to and 
mass & metallicity of its parent galaxy). The ionization field 
has little effect here, because density is the primary deter- 
minant for C II. 

C II tracing higher overdensities translates into this 
species tracing regions closer to the launch galaxy, as the 
range in d ga i in Figure |2"o1 is less than that for Civ. fi(C n) 
arises from absorbers at the same ages traced by the C iv 
absorbers (compare y-axis histograms in the center panel 
of Figure [23] and bottom left panel of I19|l . C II is tracing 
older metals that remain closer to their parent galaxies by 
virtue of either lower u w ind or more hydrodynamic slowing. 
The primary reason is the lower w W i ft d : the average t> w i n d of 
a Cn absorber is 150 kms -1 versus 250 kms -1 , while the 
Ueff/uwind ratio is not different than that for Civ; Hence ClI 
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Figure 23. Physical conditions of C II absorbers at z = 6; compare to corresponding panels for Civ in Figure Il9l and Figure |2T1 Cn is 
more dependent on density than C IV due to rapidly rising ionization fractions at higher densities. C II absorbers are more apt to trace 
metals immediately outflowing from galaxies, as well as the outflows that slow and stall closer to galaxies. Stronger C II absorbers are 
found closer to galaxies, and originate from galaxies having on average xlO less stellar mass than that for Civ absorbers. 



is arising from slower winds around less massive galaxies. 
Despite Cn tracing a different set of winds, 62% of Cn ab- 
sorbers over 50 mA are aligned with C IV in the Bubble field. 
Rare strong C n absorbers trace some of the youngest out- 
flows (10 — 30 Myrs), which are sometimes still within the 
galactic halos of z — 6 galaxies. 

In the right panel of Figure 1231 C II strength is deter- 
mined more plainly by d ga i independent of M», which is a 
result of this ion's greater dependence on density compared 
to C IV. The galaxy-absorber connection is not as important 
for Cn absorber strength. If iV(Cll) is mostly independent 
of M», the selection effects of LOSs probing a random vol- 
ume finds the average C II absorber around a smaller galaxy. 
The median M* of the £JV(C n) histogram in the bottom 
of this panel is 10 6 ' 9 M©, corresponding to a SFR of 0.03 
M0yr _1 . Cn is probing a different set of metals than Civ, 
which itself is tracing faster winds arising from more massive 
galaxies (10 7,9 Mq) enriching a larger volume. 

6.5 Evolution to z — 8 

We are already delving into the hypothetical here by con- 
sidering the physical and environmental parameters of ab- 
sorbers and galaxies not even discovered yet, so we take it 
one step further by considering what C IV absorbers trace at 
z = 8. If our suggestion that C IV absorbers are tracing met- 
als on their first journey into the IGM, metals should be ob- 
served even closer to galaxies at higher redshift. This indeed 
is the case in Figure |2"41 for z — 8 absorbers; the absorbers 
occupy slightly higher overdensity (left panel) than at z = 6 
due to winds not reaching mean overdensities yet, while ab- 
sorbers are closer to their parent galaxies as a consequence 
of having less time to travel (middle panel). These absorbers 
have u e ft more similar to « w ind signifying these outflows are 
just leaving their galaxy. In fact, the ratio w e g /f w ind for 
the absorbers between with iV(Civ) = io 13 - 5 " 14 - 5 cm" 2 for 
dgai = 10 — 20 kpc is 0.73 at z — 8 versus 0.62 at z = 6. 
Such absorbers at z — 8 trace outflowing gas more directly, 
whereas at z — 6 some absorbers at this distance trace out- 
flows that have already stalled or are returning to galaxies. 

The right panel displays our galaxy-absorber connec- 
tion plot with dotted lines corresponding to constant Bub- 
ble ionization flux, while gray lines are evolutionary tracks 



assuming v e g = 200 kms -1 and the SFR is doubled from 
Equation [4] as in DF06 at z — 8 versus z = 6. The more 
rapid assembly of galaxies at this epoch is noticeable by 
greater curvature in these evolutionary tracks; more galac- 
tic evolution likely has occurred in an absorber's associated 
galaxy despite the younger age of the Universe. The mass- 
metallicity holds at this redshift (DFO06) and appears to 
play a role in the increasing N(C iv) for more massive galax- 
ies. 



7 SUMMARY 

We explore high redshift metal-line absorbers using a 
state-of-the-art, 2 x 512 3 particle Gadget-2 cosmological 
simulation enriching the IGM with star formation-driven 
galactic outflows. Our galactic outflows follow the predicted 
relations of momentum-driven winds, where outflow velocity 
scales with galaxy velocity dispersion and the mass loading 
factor scales inversely, leading to higher mass loading from 
small galaxies. Our purpose is to reproduce IGM metal-line 
observations at z ~ 5 — 6 while providing predictions for 
future observations out to z = 8, as well as to investigate 
the physical nature of high-z metal absorbers and its 
connection to galaxies. We explore five species (Cn, Civ, 
Ol, Sill, & Siiv) using three different ionizing backgrounds: 
fully neutral, fully reionized, and a patchy bubble model 
based on the flux of the nearest galaxy. Our key results and 
predictions are as follows: 

(i) Metal lines already discovered at z > 5 are primarily 
tracing the diffuse IGM, not gas in galaxies. Some C IV lines 
already discovered could be tracing metals at overdensities 
not much higher than the mean density of the Universe. 

(ii) We can reproduce the observed trends of fi(C iv) at 
2 ~ 6 — + 4.5 with steady enrichment of the IGM filling only 
1 % of the IGM with m etals by z = 5. The recent sa mples 
of iBecker et af] [|2009h and iRvan- Weber et al l (|2009h that 
find lower f2(Clv)s suggests that the filling factor could be 
even less. We also propose that patchy ionization conditions 
where some metals are ionized by their local galaxy and oth- 
ers are not ionized may best explain the observed trends in 
both C iv and lower ionization absorbers. The determina- 



Tracing the Reionization Era with Metal Lines 31 




1 2 0.5 1 1.5 2 2.5 6 7 8 9 10 

log[p/p] log[d gal (kpc)] log[M, (M Q )] 

Figure 24. Similar to Figurc l23l but for C IV at z = 8. The color key for absorber strength from Figure [T9l applies here. See previous figure 
captions for an explanation of all contours and lines. We modify the evolutionary tracks in the right panel to reflect v e f{ = 200 kms - 1 
and double the SFR from Equation [4] reflecting galaxy evolution as found in DFO06. Absorbers still trace the location of metals well in 
this Bubble model, as shown by the overlap with metal contours in the p — T phase space of the left panel. The calculated v e ff is higher 
on average in the central panel than at z = 6. The evolutionary tracks in the right panel indicate that stronger absorbers arise from 
more massive galaxies, as at z = 6. 



tion of fl(C iv) is also highly dependent on the shape of the 
ionizing field (especially above 47.9 eV), SNe yields, the as 
parameter, and the number of N(Crv) < 10 14 cm -2 ab- 
sorbers. 

(iii) The global gas-phase metallicity increases by a fac- 
tor of 9 from z — 8 — > 5. The evolution of 0(C iv) gen- 
erally reflects this change, but the exact amount of evolu- 
tion dep ends on the form o f the io nization background. The 
uniform lHaardt fc Madaul (|200ll ) evolving background re- 
sults in an increase in Q(C iv) of x91, while the Bubble field 
shows only a x 14 increase reflecting the actual increase in 
gas-phase metals. The greater ionization intensity around 
47.9 eV in the latter background, which is a result of the 
dust attenuation law we apply along with the proximity of 
metals to galaxies, results in x8 greater f2(Clv) at Z = 8. 

(iv) C iv is the ideal tracer for metals at low overden- 
sities at z = 6, and its detection holds the best option at 
high-z to determine the volume filling factor of metals. C iv 
has the lowest global ionization correction at z ~ 5 — 6 (i.e. 
f2(C lv)/f2(C) is at its maximum). Adding a metallicity floor 
of 10 -3 Zq will have dramatic effects on the numbers of weak 
C iv absorbers, just below thresholds of currently published 
results. Integrating for several nights on a single bright z > 6 
quasar could determine if a significant portio n of the IGM 
volum e is enriched to this level. The results of iBecker et al.l 
(2009) already indicate that a significant volume of the IGM 
is unlikely to be enriched much above 10 -3 Zq, and possi- 
bly there exists too many weak C iv absorbers in our model 
without this floor. 

(v) The velocity profiles of systems are dominated by 
the peculiar velocities of outflowing gas. We reproduce the 
large Sv's (> 200 kms -1 ) observed in strong Ol and Civ 
systems. The Hubble expansion of spatially distributed met- 
als alone is far too small to create these observed system 
profiles. 

(vi) The galaxies responsible for the enrichment of the 
IGM have stellar masses in the range of lO 70 " 8 - 5 M 
and corresponding rest-frame UV magnitudes of —15.4 and 
— 18.9 (1350 A AB) when considering Civ, and a factor of 
ten less when considering C II. Most absorbers should still be 
near their originating galaxies (5-50 physical kpc for z = 6 



Civ and 2-20 kpc for Cll absorbers). We stress the need to 
identify such galaxies, which will require the likes of JWST 
and 30-meter telescope facilities, in order to study the high- 
z galaxy-absorber connection. Identifying the environments 
of strong Civ components (^ 10 14 cm -2 ) can constrain the 
ionizing source: within 30 kpc of M* = 10 70 " 90 M Q galaxies 
for the Bubble field, and at ^ 30 kpc around M» ~ 10 9 Mq 
galaxies for the HM2001 field. 

(vii) Th e abu ndance of O I absorbers observed by 
IBecker et ahl l|200rj ) in the SDSS J1148+5251 sight line is 
consistent with metals residing in a neutral IGM (as also 
argued by those authors). We cannot produce such a fre- 
quency of O I within our models with ionization extending 
above the Lyman limit. The excellent alignment with C II 
and Sill in this dataset also strongly supports this. The fact 
that most sight lines show much less Ol suggests the pos- 
sibility that reionization is still inhomogeneous over large 
scales at z > 6. It is unlikely to be an effect of inhomo- 
geneous enrichment, since smaller galaxies with relatively 
low clustering lengths are primarily responsible for enrich- 
ing the IGM at these epochs. Another possibility for broad 
Ol systems observed down to z = 5.3 are metals entrained 
in outflows, which remain shielded from the ionization back- 
ground as they travel into the IGM. 

(viii) Aligned absorber ratios are a powerful tool to de- 
termine the physical parameters of the metal-enriched IGM 
and the nature of the ionization field. Aligned absorbers 
with both components above 50 mA are ra re at z = 6, 
unless the IGM is neutral. This makes the IBecker et alj 
(2006) finding of frequent low-ionization aligned absorbers 
all the more supportive of either a neutral (probably 
recombined) portion of the IGM, or metals entrained in 
self-shielded outflows as the travel into the IGM. 

Our analysis in OD08 indicated that most metals in- 
jected into the IGM recycle back into galaxies on a typical 
timescale of 1 — 2 Gyrs. Hence at z > 5, metals should pri- 
marily be on the outward leg of their first journey into the 
IGM. This indeed turns out to be the case when we look in 
detail at the evolutionary state of absorbers. This has con- 
sequences for the galaxy- absorber connection, in that this 
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relationship is closest at high-z because metals are still in 
relatively close proximity of their parent galaxies. Absorbers 
trace gas on an outward trajectory still holding signatures 
of their galactic outflow origin, while possibly being primar- 
ily irradiated by their parent galaxy's radiation field. The 
properties of C IV absorbers depend most on their evolu- 
tionary relation to their parent galaxy (proximity, mass, & 
metallicity), while lower ionization species depend more on 
the physical state of the gas (density, primarily). This is a 
result of the different p — T phase space traced by these ions. 

Another tenet of OD08 is that metals travel a relatively 
constant physical distance (~ 100 kpc) from their parent 
galaxies, meaning early outflows have the greatest chance 
of enriching a large volume of the Universe. However this 
journey takes significant time as wind outflows travel at an 
average of ~ 100 kms -1 taking ~ 1 Gyr to get to this dis- 
tance, or more than the age of any stars at z — 6. Therefore, 
the volume of the Universe enriched to significant levels (i.e. 
0.01-0.1 Zq) grows by an order of magnitude over the rela- 
tively brief time (575 Myr) covered between z = 8 — > 5. This 
era can be called the "Age of Volume Enrichment." However 
metals do not exhibit the expansionist zeal of the Mongol 
Empire when it comes to enriching the IGM; the vast major- 
ity of the IGM (as much as 90 %) still remains unenriched in 
our simulations run to z = (|Oppenheimer fc Dave|[2009l ). 

The outlook for the study of very high-z metal lines in 
the IGM looks promising as long as sources (quasars and 
GRBs) are discovered to provide light bulbs for the near-IR 
spectrometers that will be on-line in the near future. Rela- 
tively low resolution and S/N observations can provide key 
constraints as to the filling factors of metals, the velocities 
and mass loading factors of the winds that put the met- 
als there, and of course the ionization state of the IGM. 
In contrast to metal lines between z — 5 — > 2, we expect 
significant evolution in metal-line properties between z = 
8 — > 5, a prediction that may have already been ob served by 
iBecker et all i|2009h and iRvan- Weber et al.l (|2009h . Further 
down the line is the study of the intriguing high-z galaxy- 
absorber connection, as this relation is closer than at lower 
redshift. The JWST and 20-30 meter class telescopes are 
required to study this realm as galaxies at M, < 10 s M 
with SFRs below 0.5 M©yr - , which are primarily respon- 
sible for most IGM metal lines at z ^ 6. To sum it all up, 
absorption line spectroscopy applied to the early Universe 
should yield significant insights into early galaxy formation, 
and whether such galaxies enrich the IGM in a manner sim- 
ilar to lower redshift systems, or in some exotic fashions not 
considered here. We suspect the former. 
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